AD-AI34  481  APPLICATION  OF  THERMOCHEMICAL  MODELING  TO  AIRCRAFT  i /. 

INTERIOR  POLYMERIC  MAT..(U)  JET  PROPULSION  LAB  PASADENA  » 
CA  W  DOKKO  ET  AL .  OCT  83  JPL-D-955  DOT/F AA/CT -83/ 1 6 

F/G  1/2 


UNCL'ASSI  F  I  ED 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  SUAE AU  Of  STANOASOS - ifSS " A 


Ab-AIW  V*7 


DOT/FAA/CT-83/1 6 


Application  of  Thermochemical 
Modeling  to  Aircraft  Interior 
Polymeric  Materials  II  - 
Multilayered  Seat 


W.  Dokko 
K.  Ramohalli 


Prepared  by 

Jet  Propulsion  Laboratory 
California  Institute  of  Technology 
Pasadena,  California 


October  1983 
Final  Report 


a_ 

o 

CO 


This  document  is  available  to  the  U.S.  public 
through  the  National  Technical  Information 
Service,  Springfield,  Virginia  22161. 


© 

US  Deportment  ol  Transportation 

t.  ,1  Mt|  ■  j  ,i|«|  lnl. 

rtotfw  AVKVvton  Aonwmvnmon 

Technical  Center 

Atlantic  City  Airport,  N.J.  08405 


DTIC 


SELECTE 
NOV  8  1983 


D 


83  '  l1 


07  015 


i 


i 


NOTICE 

This  document  is  disseminated  under  the  sponsorship  of 
the  Department  of  Transportation  in  the  interest  of 
information  exchange.  The  United  States  Government 
assumes  no  liability  for  the  contents  or  use  thereof. 

The  United  States  Government  does  not  endorse  products 
or  manufacturers.  Trade  or  manufacturer's  names  appear 
herein  solely  because  they  are  considered  essential  to 
the  object  of  this  report. 


1 .  Roport  No. 

DOT/FAA/CT-83/16 


4.  Titla  and  Subtitle 

Application  of  Thermochemical  Modeling  to  Aircraft 
Interior  Polymeric  Materials.il  -  Multilayered  Seat 
Cushions  _ 


7.  Author(»)  Ookko  and  Kumar  Ramohall  i 


9«  ferforvniflff  Organization  Nome  and  Addrau 

Jet  Propulsion  Laboratory 
California  Institute  of  Technology 
Pasadena.  California  91109 


12.  Sponsoring  Agoncy  Nome  and  Address 

U.S.  Department  of  Transportation 
Federal  Aviation  Administration 
Technical  Center 

Atlantic  City  Airport,  New  Jersey  08405  _ _ 


15.  Supplementary  Notes 


Technical  Report  Documentation  Page 


.  Recipient's  Catalog  No. 


5.  Roc art  Data 

October  1983 


6.  forforming  Organization  Coda 


8.  Performing  Organization  Report  No. 

JPL  D-955 


10.  Work  Unit  No.  (TRAIS) 


1 1 .  Contract  or  Grant  No. 

181-350-100 


13.  Typo  of  Roport  and  hried  Covarod 

June  1981  -  May  1982 


4.  Sponsoring  Agency  Codo 
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This  report  summarizes  the  results  from  a  twelve-month  study  of  applying  thermo¬ 
chemical  modeling  to  multilayered  polymeric  materials,  passenger  aircraft  seat 
cushions.  vThe  use  of  fire-blocking  layer (s)  between  the  foam  cushion  and  the 
covering  fabric  has  been  studied  extensively  to  minimize  fire  hazards  from  air¬ 
craft  seats.  The  objectives  of  this  work  are  to  expand  the  thermochemical  model 
for  the  multilayered  materials  and  to  experimentally  verify  theoretical  predic¬ 
tions.  First,  the  thermochemical  model  is  extended  to  any  number  of  multilayered 
materials,  by  applying  the  same  analysis  technique  used  in  the  previous  work. 

The  additional  constraints  of  temperature  and  heat  flux  continuities  at  every 
Interface  are  also  applied.  A  computer  program  is  developed  to  predict  burning 
behavior  of  seat  cushion  systems  with  and  without  a  fire-blocking  layer.  Second, 
a  series  of  tests  burning  seat  cushions  with  and  without  a  fire-blocking  layer 
are  conducted  in  a  modified  NBS  Smoke  Density  Chamber.)  The  chamber  was  supple¬ 
mented  with  a  multichannel  recorder/multiple  thermocouples  and  a  weight -measuring 
device  to  determine  temperature  profile  and  to  continuously  monitor  weight  loss, 
respectively. Vrhe  results  indicate  that  the  predicted  temperature  profiles  are 
in  very  good  agreement  with  the  experimentally  determined  ones,  and  that  the  same 
effectiveness  of  the  fire-blocking  layers  are  predicted  as  those  of  actual  weight 
loss  measurements.  It  Is,  however,  observed  that  the  formation  and  presence  of  a 
void  Inside  of  the  polyurethane  foam  seem  to  cause  the  over-prediction  of  the 
temperature  profile  and  under-prediction  of  the  weight  loss  (compared  to  the  case 
when  the  void  Is  small  or  nonexistent).  '\ 
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EXECUTIVE  SUMMARY 


The  purpose  of  this  effort  was  the  evaluation  of  the  fire  performance  of  seat 
blocking  layers  for  urethane  seat  cushions.  Because  urethane  cushions  have  been 
demonstrated  as  highly  flammable  under  aircraft  post-crash  fire  conditions,  seat 
blocking  layers  have  been  proposed  as  a  method  of  protecting  the  urethane  from 
Involvement  with  the  fire.  Many  hypothesis  have  been  proposed  to  explain  the 
manner. in  which  fire  blocking  layers  achieve  their  desired  effect,  but  such 
hypothesis  have  not  been  verified  in  any  quantitative  sense. 

A  condensed  phase  thermochemical  modeling  theory  previously  developed  for  single 
component  aircraft  materials  was  modified  to  handle  multi-layered  materials.  The 
model  is  based  on  analytical  heat  and  mass  transfer  relationships.  Required 
experimental  inputs  into  the  model  are  material  thermal  properties  typical  of  heat 
transfer  calculations.  Additionally,  the  analysis  technique  known  as  thermogravi- 
metric  analysis  is  used  to  establish  the  parameters  that  describe  the  thermal 
breakdown  of  the  plastic  material  when  exposed  to  fire. 

Parallel  to  this  analytic  work,  experimental  tests  were  conducted  on  sample  block¬ 
ing  layer  materials  in  a  modified  smoke  chamber  of  National  Bureau  of  Standards 
(NBS)  variety.  The  samples  were  exposed  to  specified  heating  rates  and  both  their 
weight  loss  and  thermal  behavior  were  measured.  This  measured  behavior  is  compared 
with  the  predictions  of  the  analytic  evaluation. 

The  findings  of  the  investigation  are  that  the  thermochemical  model  can  predict 
the  effectiveness  of  seat  blocking  layers  within  a  limited  range  of  fire  exposure 
conditions.  Additionally,  the  NBS  Smoke  Chamber  is  a  useful  small  scale  test  for 
screening  candidate  seat  blocking  layers. 


INTRODUCTION 


Aircraft  cabin  fire  safety  has  been  one  of  the  major  research  and 
development  activities  for  the  Federal  Aviation  Administration  (FAA)  for 
more  than  30  years  in  the  past  (reference  1).  The  reasons  for  such  a 
long-term  commitment  are  the  fatalities  observed,  unique  environment  in  an 
aircraft  (such  as  small  space,  high  density  of  people,  limited  access  for 
egress,  etc.),  and  the  complexity  in  fire  phenomenon  itself. 

In  general,  fire  potential  in  a  passenger  aircraft  can  be  due  to  jet  fuel 
and  cabin  interior  materials.  However,  a  fire  involving  only  the  interior 
materials  can  occur  in-flight  or  on  the  ground.  In  the  case  of  a  postcrash 
fire,  a  jet  fuel  fire  can  ignite  interior  materials,  and  the  hazards  arise 
from  interaction  of  the  fuel  fire  with  interior  materials  (reference  2). 

As  planned  earlier  (reference  4),  this  work  is  on  the  prediction  of 
material  burning  phenomena  through  condensed  phase  thermochemical  modeling, 
and  does  not  involve  gas  dynamic  modeling.  These  two  approaches  may  well 
be  bridged  in  the  future  (reference  3)  with  sufficient  progress. 

Among  the  cabin  interior  materials  such  as  carpet,  seats,  window  screens, 
sidewall  panels,  ceiling  and  partitions,  one  of  the  most  flammable 
materials,  seat  cushion,  is  the  subject  of  the  current  work.  Previous 
work,  last  year,  (reference  4)  was  on  wool  carpet  and  polyurethane  foam 
treated  as  single  layer  materials. 

In  this  work,  the  seat  is  treated  as  a  system  of  multilayered  polymeric 
materials,  consisting  of  a  seat  cover  fabric,  polyurethane  foam  cushion, 
and  a  fire  blocking  layer. 

In  addition  to  analytically  predicting  the  burning  behavior  of  multilayered 
systems  as  a  function  of  heat  flux  and  layers'  thicknesses,  burning  tests 
are  conducted  in  a  modified  NBS  density  chamber  to  verify  the  temperature 
profile  and  weight  losses  predicted  by  the  model. 


2.  THEORETICAL 


The  thermochemical  model  used  in  this  work  is  an  extension  of  the  model 
developed  earlier  for  a  single-layered  material.  As  a  matter  of  fact,  the 
approach  used  in  predicting  the  burning  rate  of  wool  carpet  with  a  char 
layer  on  its  top  (reference  4)  is  a  particular  form  of  such  an  extension. 
Another  earlier  work  (reference  5)  applied  to  the  study  of  burning  rate  of 
composite  materials  is  also  referenced  for  the  generalized  extension. 

(a)  The  Applicability  of  the  Steady-State  Model 

There  are  two  aspects  of  this  time-independent  model  which  shot  '  be 
recognized.  One  is  the  coordinate  transformation  (a  la  Spaldi 
(reference  11))  that  enables  us  to  treat  the  regressing  surfac  •  if 
it  were  stationary  and  the  other  is  the  more  fundamental  arguir  on 
the  relative  time  scales. 

(1)  The  assumption  of  time-independent  degradation  enables  one  to  see 
that  d/dt  can  be  written  as  r  •  d/dx.  This  transformation,  the 
authors  believe,  was  first  introduced  in  combustion  theory  by 
Brian  Spalding  (reference  11)  in  his  treatment  of  laminar  flame 
propagation  in  premixed  gases.  The  point  is  that  the  “transient" 
propagation  can  be  viewed  as  steady-state  in  another  coordinate 
frame. 

(2)  In  any  case,  the  assumption  8/9  t  =  0  implies  that  the  charac¬ 
teristic  heat  transfer  time  is  small  compared  to  the  characteris¬ 
tic  "flow"  or  thermal  wave  propagation  time.  The  time  needed  for 
the  establishment  of  the  fully  developed  temperature  profile  in 
the  solid  is,  of  course,  infinity  from  the  first  order  nature  of 
the  equation's;  recall 

(Tx  -  T0)/(TXS  -  T0)  =  1  -  exp  (-ot/x?) 

where  Tx  is  the  temperature  at  any  depth  x  (from  the  surface), 

T0  is  the  initial  temperature,  Txs  is  the  steady-state 
temperature  at  depth  x,  is  the  thermal  diffusivity 
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(  s  k/  pc),  t  is  the  time.  Thus  the  time-independence 
assumption  in  our  treatment  relies  on  a  small  value  for  the 
characteristic  time  l^/a  for  its  validity.  In  the  case  on 
hand,  the  following  approximate  numbers  lend  credence  to  this 
approach.  Typically  a  -  10“^  cm^/sec.  The  "characteris¬ 
tic"  thermal  depth,  i.e.,  the  distance  from  the  surface  to  reach 
1/e  of  the  full  temperature  difference  is  0.2  cm  (for  example 
from  figure  11,  first  plot)  at  120  seconds  (2  min.);  or 

t  -iL  =3.0 
x^ 


leading  to 

1  -  exp  (-  at/x 2)  >  0. 9 

Thus  the  steady-state  assumption  is  valid  in  our  case  of  heating 
rates. 


(b)  Mathematical  Formulation 

As  in  the  previous  work,  the  geometry  under  consideration  is  one¬ 
dimensional,  and  the  steady-state  condition  is  assumed  throughout  the 
analysis.  The  following  mathematical  analyses  for  a  multilayered 
system  are  divided  into  three  parts:  the  first  part  is  for  the  case 
where  a  single-layered  material  is  so  thick  that  the  rear  surface 
(cold  side)  temperature  is  undisturbed;  the  second  part  is  the  same  as 
above  except  that  the  rear  surface  temperature  has  been  raised 
substantially,  and  the  third  case  is  the  extension  for  a  layer  which 
is  on  top  of  the  above  two  layers  and  has  as  many  colder  layers  under¬ 
neath  it  as  determined  by  the  configuration. 
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In  all  three  cases,  the  governing  equation  for  each  layer  is  the  same, 
that  is, 

k  flil  +  per  £L  =  DpNB  exp  (-E/RT)  (1 ) 

dx2  dx 

The  right-hand  side  of  this  equation  is  the  heat  sink  term 
corresponding  to  the  first-order  degradation  reaction  described  by 

1  dN 

=  B  exp  (-E/RT)  (2) 

The  symbol  N  represents  the  number  of  remaining  bonds  per  unit  mass, 
normalized  with  respect  to  that  of  fresh,  unburned  material.  The 
symbol  D,  the  heat  of  degradation,  is  "positive"  when  the  degradation 
reaction  is  endothermic  and  "negative"  when  exothermic. 

The  boundary  conditions,  however,  are  different  for  each  case. 

(1)  First  case,  i.e.,  T  =  T0  at  x  =  °o 


Figure  1.  Geometry  considered  in  the  first  case  where  T  *  T0  at  x  «  oo 
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When  the  rear  surface  (cold  side)  temperature  is  undisturbed,  the 
boundary  conditions  for  equations  (2)  and  (3)  are,  respectively. 


T 

*  Ts 

at 

O 

II 

X 

(3) 

T 

=  To 

at 

x  :  •» 

N 

=  1  - 

o 

II 

X 

(4) 

FS 

N 

=  1 

at 

X  =  oo 

The  symbol  FS  stands  for  fragment  size,  an  average  number  of  monomer 
units  of  degraded  polymer. 

Assuming  constant  values  of  material  properties,  the  solution  for  r, 
surface  regression  rate,  has  been  obtained  (reference  6)  by  singular 
perturbation  methods  with  matching  the  solution  for  the  inner 
(surface)  and  outer  (deep  solid)  regions.  For  the  steady  state 
system,  r  is  given  as: 


r  = 

And  the  integration  of  equation  (1)  will  show  that  the  conductive  heat 
flux  at  the  hot  surface  is 


(k/pc)  B  exp  (-E/RTJ 


am 


c(Ts-T0)FSj 


(5) 


qs  =  pr  [c  (Ts-T0)  +  d/fs3  (6) 

if  there  is  zero  heat  flux  at  the  rear  surface.  Thus,  all  the  thermal 
energy  flowing  into  the  layer  can  be  said  to  be  totally  consumed 
within  the  layer. 
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(2)  Second  Case,  i.e.,  T  =  Tc  (>  T0)  at  x  >  0 


Figure  2.  Geometry  considered  In  the  second  case  where  T  ■  Tc  at 
x  *  x^ 


This  Is  the  case  when  the  front  (hot)  and  rear  surface  tempera¬ 
tures  of  a  layer  are  substantially  higher  than  the  unperturbed 
condition  of  T0.  The  governing  equations  are  the  same  as 
above,  i.e.,  equations  (1),  (2),  and  the  boundary  conditions  are 


T  *  Ts  at  x  =  0 

T  =  Tc  at  x  =  xL 


N  *  1  -  -7  at  x  =  0 

FSs 

N  =  1  -  77  at  x  =  x. 


(7) 

(8) 


Assuming  that  Tc  is  close  to  Ts,  that  is 


111 

VT 


c 


0 


<  0.1 


(9) 


the  solution  for  r  is  given  by 


r 


(k/pc)  B  exp  ( -E/RT s ) 

-  go  m 


(10) 
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The  details  of  this  derivation  Is  shown  In  appendix  A. 

The  Integration  of  equation  (1)  In  the  region  between  the  front 
and  rear  surfaces  will  give  the  net  heat  flux,  that  Is,  the 
Incoming  heat  flux  at  the  front  surface  minus  the  one  at  the  rear 
surface. 

qnet  =  pr  <VTC>  +  d{,fs7  ~ 

In  other  words,  the  above  equation  determines  the  net  heat 
consumed  within  the  layer  volume,  as  sensible  heat  and  heat  of 
degradation. 

(3)  Third  Case  of  a  multilayered  system 

In  the  above  two  cases,  the  analyses  are  for  a  homogeneous  single 
layer  for  which  thermal,  physical  and  chemical  properties  are 
constant  and  uniform  in  every  part  of  the  layer.  If  a  system 
consists  of  more  than  one  material,  like  the  example  of  a  cover 
fabric  -  blocking  layer  -  foam,  those  properties  are  widely 
different  and  hence  will  be  treated  as  a  multilayered  system. 

Strictly  speaking,  even  a  single  component  material  should  be 
treated  as  a  multilayered  system,  if  the  temperature  drop  from 
the  front  to  the  rear  surface  is  so  large  that  the  thermal, 
physical  and  chemical  properties  can  not  be  considered  constant 
and  uniform  within  the  layer. 

In  other  words,  when  a  material  of  certain  thickness  is  burning, 
it  is  a  multilayered  system  from  the  view  point  of  thermochemical 
behavior.  As  shown  later,  B  and  E  of  any  single  material 
employed  in  this  study  do  not  remain  constant  over  the  tempera¬ 
ture  range  of  interest  but  vary  considerably. 
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For  this  reason,  each  component  material  of  a  multilayered  seat 
cushion  system  Is  In  Itself  considered  multilayered,  divided  by 
temperatures  at  which  any  of  the  thermal  and/or  thermochemical 
properties  changes  substantially.  Hence,  the  total  number  of 
layers,  l.e.,  thermochemical  layers,  of  a  seat  cushion  system 
become  equal  to  the  sum  of  the  thermochemical  layers  of  each 
component  material. 

For  the  ease  of  analysis,  first  suppose  that  a  system  of  three 
layers  Is  undergoing  pyrolysis  reaction,  as  depicted  In  figure  3. 
The  same  analysis  can  be  extended  to  any  system  with  more  than 
three  layers  with  the  change  in  subscripts. 


Figure  3.  Geometry  of  the  third  case,  a  three-layered  system 

The  conductive  heat  flux  at  the  top  surface  (x*0)  should  be  the 
sum  of  all  the  energy,  i.e.,  the  sensible  heat  and  the  heat  of 
degradation,  required  for  the  pyrolysis  of  all  three  layers  In 
this  1-0  approximation.  The  energy  required  for  each  Individual 
layer  Is  shown  In  equation  (6)  or  equation  (11).  Then,  the 
Incoming  heat  flux  at  the  top  surface  of  the  first,  second,  and 
the  third  layer  is,  respectively. 
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The  above  two  conditions  are  to  satisfy  the  continuity  of 
temperature  and  heat  flux  at  interfaces. 

In  the  expression  of  equations  (12),  (13),  or  (14),  each  r  is 
implicitly  defined.  In  other  words,  the  heat  fluxes  at  the  upper 
and  lower  boundaries  of  each  layer  are  used  in  obtaining  the  r  of 
that  layer  (as  shown  in  appendix  A),  and  are  in  turn  used  to 
define  the  heat  fluxes  [as  shown  in  equations  (12),  (13),  and 
(14)].  This  poses  no  problem  in  actual  calculations  since  any 
iterative  method  leads  to  fairly  rapidly  converging  answers.  The 
significance  of  this  fact  is,  however,  that  the  calculation  of 
r's  should  proceed  layer  by  layer,  from  the  rear  (coldest)  to  the 
front  (hottest),  or  from  the  front  to  the  rear. 

(c)  Computer  Program 

Based  on  the  theory  discussed  above,  two  computer  programs  are 
developed  and  listed  in  appendix  B.  The  first  program,  (Computer 
Program  A)  takes  the  interfacial  temperatures  as  the  input  data,  and 
the  second  (Computer  Program  B)  takes  the  temperature  and  the  conduc¬ 
tive  heat  flux  at  the  front  surface  as  its  input.  The  choice  among 
these  two  programs  is  dictated  firstly  by  the  available  boundary  con¬ 
ditions,  and  the  other  may  be  used  as  a  complementary  to  obtain 
additional  information. 

(d)  The  Overall  Logic  and  Characteristics  of  the  Model 

The  model  was  originally  developed  to  account  for  the  host  of 
inconsistencies  in  polymer  degradation  and  burning  rate  data.  The 
most  important  distinguishing  characteristics  are  the  recognition  that 
the  mean  molecular  weight  of  the  vaporizing  molecules  need  not  be 
equal  to  the  monomer  molecular  weight.  Instead,  a  vapor  pressure 
equilibrium  criterion  was  used  to  unambiguously  specify  the  mean 
molecular  weight  at  the  surface.  This  procedure  not  only  removed  an 
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annoying  arbitrainess  at  the  surface,  but  also  checked  well,  in  terms 
of  results,  with  data  from  various  classes  of  polymer  combustion. 

Under  the  assumption  of  first  order  Arrhenius  kinetics  for  the  degra¬ 
dation  of  the  subsurface  polymer,  the  model,  in  its  present  state, 
cannot  handle  cases  vrfiere  extensive  degradation  beyond  the  monomer 
stage  takes  place  before  vaporization.  This  is  not  as  severe  a 
limitation  as  it  may  seem.  It  has  been  documented  (e.g.,  see  Stanley 
Martin,  X  Symposium  (International)  on  Combustion)  that  many  of  the 
the  degradation  reactions  actually  take  place  after  the  major  precur¬ 
sors  (i.e.,  products  of  partial  degradation)  have  left  the  surface. 

The  model  also  assumes  that  the  thermal  wave  front  moves  at  a  uniform 
speed  in  the  solid.  This  thermal  wave  can  be  associated  with  the  mean 
surface  (and  hence  the  wave  speed  is  indeed  the  burning  rate  or 
"regression  rate")  only  when  a.  5  of  the  solid  (condensed  phase)  gets 
transformed  into  vapor.  This  is  indeed  the  case  in  the  burning  of 
simple  plastics  (e.g.,  plexiglas,  polyethylene,  polyurethane  .  .  .  ). 
In  the  case  on  hand,  substantial  charring  occurs  indicating  that  the 
familiar  "regression  rate"  needs  careful  interpretation.  For  example, 
it  is  a  familiar  fact  that  in  wood  burning  the  reactive  portion 
(cellulose  and  hemicell ulose)  leaves  the  solid  framework  of  lignin 
char.  In  such  cases  of  charring  solids,  the  regression  rate  is  more 
appropriately  associated  with  the  velocity  of  the  thermal  wave  in  the 
wake  of  which  the  reactive  portion  gets  vaporized  and  leaves  the  char. 
Thus,  it  may  appear  that  the  surface  is  not  really  regressing  in  the 
physical  sense,  although  the  reactive  portion  is.  Up  to  this  point, 
the  interpretation  of  the  regression  rate  of  this  model  is  clear.  The 
interpretation  gets  more  complex  as  multilayered  materials  are 
concerned.  Hence,  clearly,  the  regression  rate  of  the  XT  layer,  if 
interpreted  literally,  will  give  rise  to  a  "void"  or  separation 
between  the  bottom  of  the  I  layer  and  the  initial  position  of  the  XL 
layer  top  surface.  However,  the  complexity  is  removed  when  one 
recognizes  that  the  regression  meant  here  is  really  the  movement  of 
the  thermal  wave  in  the  solid  layers  vrtnle  the  char  framework  (end 
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product  of  the  degradation)  is  stationary  in  the  laboratory  frame  of 
reference.  The  appropriate  way  of  handling  this  situation  mathemati¬ 
cally  is  to  keep  the  density  of  the  material  P  as  a  variable.  For 
example,  P  final  would  not  equal  2ero  but  would  equal  the 
char  density  which  itself  is  equal  to  the  (P  initial  — 

Preactive*) 

With  this  interpretation,  the  regression  rate  r  is  really  always 
associated  with  the  mass  flux  p  r  where  P  is  not  the  density  of  the 
solid  but  is  only  the  density  of  the  reactive  portion  of  the  solid. 
However,  in  the  evaluation  of  the  properties  such  as  the  thermal 
diffusivity  of  the  reacting  solid,  the  full  density  has  to  be  used. 
Also,  the  heat  transfer  through  the  char  is  characterized  by  a  thermal 
conductivity  coefficient  of  the  char,  and  the  volatiles  convection 
(flow)  through  the  char.  This  has  been  highly  simplified  in  the 
present  analysis  which  considers  the  temperature  to  be  constant  (and 
equal  to  the  temperature  of  the  bottom  surface  of  the  immediate  layer 
above).  Clearly,  more  work  is  needed  to  mathematically  incorporate 
into  the  model,  the  realistic  char  formation,  flow  and  heat  transfer. 

The  "void"  that  is  mentioned  in  the  experimental  portion  is  not  the 
same  on  this  void  or  separation.  Experimentally  it  is  found  that  the 
mechanical  movement  of  the  volatiles  frequently  causes  a  void  in  the 
assembly  and  physical  separation  of  layers.  This  is  a  very  complex 
process  mathematically. 

Considering  all  of  these  complexities,  the  agreement  between  the 
theory  and  experiments  in  this  study  is  thought  to  be  encouraging. 

3.  EXPERIMENTS 

The  experiments  of  this  work  consist  of  two  parts;  one  is  to  obtain  the 
kinetics  constants  using  thermogravimetry  (TG)  and  the  other  is,  through 
actual  burning  tests,  to  determine  the  temperature  profile  established 
within  a  seat  cushion  assembly  and  the  weight  loss  during  pyrolysis.  The 
experimental  results  are  compared  with  the  model  predictions. 
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(a)  Thermogravimetry 


This  experiment  heats  up  a  small  sample  suspended  inside  a  furnace  and 
records  the  weight  change  as  a  function  of  temperature.  The  kinetics 
constants,  B  and  E,  are  deduced  then  from  an  Arrhenius  plot,  a  plot  of 
log  j(l/*/)  •  -(dw/dt)}  versus  reciprocal  of  absolute  temperature 
O/j)* 

It  was  known  from  previous  work  (reference  4)  that  the  ambient  gas 
composition,  more  precisely  the  oxygen  concentration,  affects  the 
kinetics  constants.  This  is  a  very  important  point.  It  has  long  been 
recognized  that  even  minor  ("■'1%)  concentrations  of  certain  oxidative 
species  can  significantly  alter  (by  almost  an  order-of-magnitude)  the 
degradation  rates  of  polymers.  And  yet  no  detailed  study  seems  to  be 
available  in  this  area.  One  study  that  addressed  this  question 
specifically  stopped  short  of  actually  demonstrating  the  effects  but 
had  to  make  valid  (but  indirect)  deductive  arguments  to  point  the 
importance.  In  any  combustion  situation,  especially  with  flow  of 
gases  over  the  surface  (a  very  common  aircraft  fire  scenario),  the 
oxygen  concentration  £t  the  burning  surface  appears  to  be  around 
0. 1%-1%  (Wooldridge  and  Muzzy  (reference  8),  Kulgein  (reference  9), 
Fennimore  and  Jones  (reference  10).  This  measured  non-zero 
concentration  is  significant,  because  earlier  JPL  work  under  FAA 
sponsorship  (reference  4)  indicated  that  actual  aircraft  interior 
materials  exhibited  a  strong  dependence  of  the  kinetics  constants  of 
degradation  on  small  concentrations  of  oxygen  in  a  stream  of  oxygen  in 
an  inert.  When  it  is  realized  that  turbulent  transport  offers  a 
mechanism  for  the  availability  of  small  concentrations  of  oxygen  at 
the  burning  surface,  it  is  easy  to  see  that  these  oxygen  concentration 
effects  could  be  important  in  predicting  full-scale  burning  behavior. 
Full-scale  diffusion  flames  would  be  large  enough  to  be  fully 
turbulent. 

Nevertheless,  all  the  present  TG  experiments  are  performed  in  pure 
nitrogen  except  for  seat  cover  fabric,  due  to  the  limited  scope  of  the 
work. 
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The  T6  diagrams  of  the  cover  fabric,  fire  blocking  layers  (Vonar®, 
LS-200,  and  Preox®  ),  and  foams  (polyurethane  and  polyimide)  are 
shown  in  figures  4,  5,  and  6,  respectively.  The  Arrhenius  plots  of 
figures  7,  8,  and  9  show  that  the  pyrolysis  of  each  of  these  materials 
cannot  be  described  as  a  single  first-order  reaction.  Thus,  the 
entire  reaction  temperature  range  is  subdivided  into  multiple  segments 
in  each  of  which  the  reaction  rate  is  reasonably  accurately  described 
by  the  first-order  Arrhenius  expression.  Figure  7  shows  how  it  is 
done  for  the  cover  fabric,  as  an  example,  with  six  subdivided 
segments.  Each  segment  comprises  one  thermochemical  layer  as 
discussed  in  the  previous  section. 

Kinetics  constants  obtained  along  with  other  thermal  and  physical  data 
of  the  materials  used  in  this  study  are  shown  in  Table  I.  These  are 
used  as  input  data  for  the  computer  program  shown  in  Appendix  B. 

(b)  Burning  Tests 

In  order  to  measure  the  temperature  profile  and  weight  loss  of 
samples,  a  modified  NBS  Smoke  Density  Chamber  is  used.  The  radiative 
heat  flux  was  provided  by  a  high  heat  flux  furnace  (Mellen  furnace 
Model  No.  10,  with  the  maximum  heat  flux  of  12  W/cm2  or  10.6 
BTU/ft2  sec).  The  originally  equipped  furnace  for  the  NBS  chamber 
could  provide  2.5  W/cm2.  For  the  temperature  distribution  measure¬ 
ment,  nine  thermocouples  (Pt  vs.  Pt-10%  Rh)  were  connected  to  a 
multichannel  recorder  (Leeds  and  Northrup's  Speedomax  Model  251,  12 
channel,  with  resolution  time  of  1  second  for  each  reading).  The 
weight  measuring  device  was  a  transducer-type  cantilever  beam, 
connected  to  an  amplifier  and  then  to  a  strip-chart  recorder. 

The  furnace  is  heated  up  to,  and  maintained  at,  850°C  during  the 
tests.  After  the  furnace  has  reached  a  steady  state,  a  radiation 
shield  is  removed  and  the  seat  cushion  is  exposed  directly  to  the 
furnace.  The  radiative  heat  flux  is  around  ~1.6  Btu/ft2  sec. 
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FIGURE  4.  THERMOGRAVIMETRIC  DIAGRAM  OF  SEAT  CUSHION  COVER  FABRIC 


FIGURE  5.  THERMOGRAVIMETRIC  DIAGRAMS  OF  THREE  FIRE  BLOCKING  LAYERS,  LS-?00,  VONAR,  AND  PREOX 
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FIGURE  6.  THERMOGRAVIMETRIC  DIAGRAMS  OF  TWO  FOAMS,  POLYURETHANE  FOAM  AND  POLYIMIDE  FOAM. 


FIGURE  8.  ARRHENIUS  PLOTS  OF  TWO  FIRE  BLOCKING  LAYERS,  VONAR  AND  PREOX, 
DEDUCED  FROM  TG  SHOWN  IN  FIGURE  5. 


Table  1(a).  Fundamental  Properties  of  Seat  Cushion  Component  Material 


c 

k* 

*L 

Density 

Thermal 

9 

Specific  Heat 
cal 

Conductivity 

cal 

Thickness 

cm  ^ 

q.at 

cm.  sec.°C 

cm 

Seat  Cover  Fabric 

0.4 

0.3 

0.0001 

0.1 

Fire  Blocking  Layer 

LS-200 

0.12 

0.3 

0.0002 

1.2 

Vonar® 

0.146 

0.345 

0.0002 

0.8 

Preox® 

0.62 

0.3 

0.00034 

0.1 

Foam 

Polyurethane 

0.03 

0.4 

0.0001CK0.00034 

10.6 

Polyimide 

0.023 

0.2 

0.0001 

(b)  Kinetics  Constants  and  Heat  of  Degradation 
of  Subdivided  Thermochemical  Layers 


T 

Tempera¬ 

ture 

Range 

(K) 

B 

Pre  Exponential 
Factor 

(1/sec) 

E 

Activation 

Energy 

(cal/mole) 

D 

Heat  of 
Degrada¬ 
tion 

(cal/g) 

848-898 

2.13  x  104 

24800 

-2.9 

Seat 

813-848 

5.06  x  10"8 

-20300 

-2.9 

Cover 

763-813 

2.24  x  1013 

56500 

-2.9 

Fabric 

689-763 

2.51  x  10-1 

-2.3 

602-689 

7.30  x  10"6 

-6500 

-16. 

523-602 

1.57  x  102 

13700 

-16.3 

*  Source:  NASA  Ames  Research  Center  —  Handbook  of  Chemistry  ond Physics 
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(b)  Kinetics  Constants  and  Heat  of  Degradation 
of  Subdivided  Thermochemlcal  Layers  (Cont'd) 


T 

B 

E 

D 

Tempera- 

Heat  of 

ture 

Pre  Exponential 

Acti vat  ion 

Degrada- 

Range 

Factor 

Energy 

tion 

(K) 

(1/sec) 

(cal /mole) 

(cal/g) 

877-923 

1  1.17  x  10-5 

-3857 

0.5 

767-877 

1  7.17  x  10"11 

-247/0 

9.6 

723-767 

5.17  x  101 

16818 

-8.4 

698-723 

1.04  x  10"}. 

-2020 

-19.1 

LS-200 

673-698 

8.28  x  10%4 

-31083 

-16.7 

648-673 

6.56  x  10;3 

2475 

-8.4 

623-648 

2.46  x  106 

27897 

-20.3 

598-623 

1.61  x  10;17 

-38190 

7.9 

573-598 

2.11  x  107 

27800 

71.9 

548-573 

2.01  x  10:1 

6772 

24.9 

473-548 

7.83  x  103 

18279 

54.9 

881-923 

3.37  x  10‘4 

1043 

1.0 

766-881 

1.66  x  10~7 

-12285 

9.6 

709-766 

8.94  x  10"f, 

7800 

-23.9 

661-709 

2.37  x  10"lz 

-26532 

-28.7 

Vonar^ 

623-661 

2.75  x  104 

22078 

-28.7 

591-623 

3.94  x  10"13 

-25933 

50.9 

473-591 

2.33  x  104 

19435 

107.5 

873-908 

5.42  x  10;6 

-13061 

0.5 

848-873 

8.55  x  1020 

91578 

0.5 

Preox ^ 

748-848 

1.55  x  10% 

18860 

0.5 

683-748 

9.85  x  10'3 

4500 

-13.9 

623-683 

1.43  x  10-5 

-4365 

-13.9 

573-623 

5.87  x  10;3 

3085 

-13.9 

498-573 

5.94  x  10° 

10963 

0.5 

Polyurethane 

644-688 

8.44  x  1019 

65760 

3.1 

Foam 

598-644 

6.17  x  101 

12335 

11.9 

573-598 

1.77  x  10;2 

2646  1 

2.5 

544-573 

2.60  x  105 

21435  | 

2.5 

448-544 

6.00  x  10° 

9900 

0.5 

1006-1073 

3.06  x  10-5  | 

-4763 

0.5 

973-1006 

3.12  x  10- J4  1 

-46170 

0.5 

Polyimide 

915-973 

1.80  x  10-7  1 

-16070  , 

0.5 

Foam 

854-915 

1.33  x  10% 

16860 

-7.3 

796-854 

6.11  x  10-5 

-7890 

-7.3 

711-796 

1.07  x  107 

36705 

673-711 

7.63  x  IQ"3 

6938 

-0.3 
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(c)  A  Note  on  the  Heating  Rates 


The  actual  heat  fluxes  encountered  in  aircraft  crash  situations  can 
approach  10  W/cm^  —  a  fact  that  led  investigators  to  recognize  that 
the  NBS  smoke-density  chamber,  with  its  original  2.5  W/m^  furnace, 
may  have  to  be  modified  for  realistic  simulations.  In  the  present 
effort,  the  heat  fluxes  meant  for  the  FAA  samples  were  planned  to  be 
gradually  increased  from  a  low  value  (1  W/cm^)  to  the  full  value  of 
6-10  W/cm^.  It  was  the  plan  to  verify  the  model  at  low  heat  fluxes 
first  in  a  logical  sequence  of  increasing  heat  fluxes.  Thus,  at  the 
time  of  the  present  reporting  the  lower  heat  fluxes  have  been 
investigated.  The  capability  exists  for  the  higher  heat  fluxes. 
However,  the  actual  tests  await  future  research  support. 

Since  the  initial  setup  of  equipment  and  trial  runs,  several  modifica¬ 
tions  and  improvements  have  been  made  on  the  sample  preparation  and 
run  procedure,  in  order  to  obtain  meaningful  and  reproducible  test 
results.  These  are  listed  in  appendix  C  for  future  reference. 

A  pair  of  thermocouples  are  placed  at  every  interface,  including  the 
front  surface.  And  additional  five  are  placed  within  the  polyurethane 
foam,  0.5  inches  apart  from  each  other  along  the  axis.  A  test  is 
always  run  in  duplicate;  one  to  measure  the  weight  loss  and  the  other 
for  the  temperature  profile.  This  is  because  the  thermocouples  have 
such  stiffness  that  their  use  would  not  allow  free  movement  of  the 
cantilever  beam  and  the  accurate  measurement  of  weight  change. 

RESULTS  AND  DISCUSSION 

(a)  Burning  Test  Results 

The  pictures  of  the  burned  samples  and  the  equipment  used  in  this  work 
are  shown  in  appendix  D  for  visual  examination.  The  results  of  the 
burning  tests  in  the  modified  NBS  Smoke  Density  Chamber  are 
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shown  in  figure  10.  These  are  the  temperature  changes  obtained  from 
the  multiple  thermocouple  readings  placed  at  different  positions 
within  the  seat  cushion  assemblies  and  the  weight  change  with  time, 
for  the  runs  with  and  without  the  fire  blocking  layers. 

The  weight  change  data  show  that  the  sample  without  a  fire  blocking 
layer  suffered  the  largest  mass  loss,  and  that  the  effectiveness  of  a 
fire  blocking  layer,  if  judged  from  the  mass  loss  data,  can  be  rated 
in  increasing  effectiveness  on  PreoxdD,  VonarCD,  and  LS-200.  The  same 
order  is  observed  when  the  burnt  samples  were  examined  visually. 


The  thermocouple  reading  at  the  sample  surface  facing  the  furnace, 
however,  should  be  corrected  for  the  caused  by  radiation  and  convec¬ 
tion.  The  energy  balance  for  the  thermocouple  bead  can  be  shown  as 


*Yad  ^furnace  ”  ^bead^ 


h  (T.  .  -  T  ...  ) 

conv  '  bead  gas  film' 


where  hracj  and  hconv  are  the  radiative  and  convective  heat  transfer 
coefficient,  respectively,  and  Tfurnace,  Tbead,  Tgas  film  are  the 
temperatures  of  the  furnace,  thermocouple  bead,  and  film  stagnant 
gas  film  surrounding  the  bead,  respectively.  The  estimated  error 
(^actual  ~  Tbead)  is  ’n  t,1e  range  of  25°C  to  35°C  under  the  current 
experimental  conditions.  Further  detailed  error  analysis  is  not 
attempted  here  due  to  the  lack  of  precise  information  on  emittance, 
sample  surface  condition,  and  thermal  contact  between  the  bead  and 
fabric  fibers,  which  are  necessary  to  estimate  the  heat  transfer 
coefficients  with  sufficient  precision. 


Other  than  this,  the  thermocouple  data  lead  to  a  few  important 
observations. 

(1)  The  first  one  is  that  the  front  surface  temperature  measured  by 
the  thermocouple  No.  1  reaches  a  steady  state  in  about  three 
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FIGURE  10. 

FIGURE  10.  WEIGHT  AND  TEMPERATURE  CHANGE  DATA  MEASURED  DURING  THE  BURN 
TESTS  IN  A  MODIFIED  NBS  S^OKE  CHAMBER  WITH  AND  WITHOUT  A 
FIRE  BLOCKING  LAYER  (2  of  2  pages) 
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minutes,  and  then  registers  essentially  the  same  temperature 
thereafter.  When  a  fire  blocking  layer  of  Vona  r®  or  LS-200  is 
used,  however,  a  very  slight  increase  of  the  front  surface 
temperature  is  observed.  The  temperature  just  behind  the  cover 
fabric,  measured  by  the  thermocouple  No.  2,  trails  that  of  the 
front  surface  by  about  15°C. 

(2)  There  is  a  further  temperature  drop  across  a  fire  blocking  layer. 
This  temperature  drop  is  indicated  by  the  vertical  distance 
between  the  thermocouples  No.  2  and  No.  3  readings,  as  shown  in 
figure  10.  In  the  figure,  LS-200  allows  the  largest  temperature 
drop,  Vonar® medium,  and  Preox® the  least.  Incidentally 

this  is  exactly  the  opposite  order  as  the  one  observed  in  total 
mass  loss. 

(3)  The  temperature-rise  histories  within  the  polyurethane  foam, 
measured  by  thermocouples  No.  4,  5,  6,  and  7,  hardly  represent 
steady-states.  Most  of  the  measured  temperatures  simply  keep 
rising  steadily,  as  shown  in  figure  10. 

In  some  cases,  however,  a  plateau  seems  to  appear  after  such 
monotonic  increase  in  temperature,  implying  a  steady -state  condi¬ 
tion  attained  within  the  foam.  This  behavior  is  believed  to  be 
caused  by  the  formation  and  existence  of  a  void  or  dome  filled 
with  relatively  hot  pyrolysis  gas  products  in  a  convective  flow 
motion.  This  explains  why  such  a  plateau  shows  up  and  why  then 
these  plateaus  are  close  to  each  other  -  in  other  words,  the 
temperature  gradient  is  small.  It  is  particularly  serious  when 
no  fire  blocking  layer  is  used,  because  under  the  condition  of 
high  heat  flux,  the  polyurethane  foam  liquefies  and  drips  before 
complete  burnout  occurs.  This  eventually  leaves  an  ever-growing 
void  filled  with  the  hot  pyrolysis  gas  product.  The  hot  gas  of 
the  void  may  flow  freely  and  contribute  to  the  rapid  temperature 
rise  at  the  rear  surface. 
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(4)  The  use  of  a  fire  blocking  layer  helps  maintain  lower 

temperatures  at  the  rear  (colder)  side  essentially  by  reducing 
the  heat  transfer.  This  can  be  seen  easily  by  comparing  the 
temperatures  measured  by  the  thermocouples  with  the  same  number. 
For  example,  after  6  minutes  of  exposure,  the  thermocouple  No.  5 
reaches  321  K  with  LS-200,  323  K  with  Vonar®,  and  463  K  with 
Preox®,  while  it  reaches  546  K  without  a  fire  blocking  layer. 

There  is,  however,  a  compensating  effect  near  the  front  (hot) 
surface.  That  is,  the  front  attains  higher  temperatures  with  the 
fire  blocking  layer  than  would  be  observed  without  the  fire 
blocking  layer.  Figure  10  shows  the  front  surface  temperature  is 
701  K  with  LS-200,  695  K  with  Vonar®,  691  K  with  Preox®, 
and  643  K  without  one,  respectively,  after  10  minutes  exposure. 
This  leads  to  an  observation  that  a  fire  blocking  layer  functions 
as  if  it  is  a  reflector  for  heat  flux. 

This  results  in  lowering  the  rear  side  temperature  while 
maintaining  the  higher  front  side  temperature. 

(b)  Predictions  by  Thermochemical  Model 

The  temperature  profiles  obtained  by  the  thermochemical  model,  in 
comparison  with  the  experimental  data,  are  shown  in  figure  11.  The 
profile  for  each  configuration  (with  and  without  a  fire  blocking 
layer)  is  calculated  every  2  minutes  using  the  Computer  Program  B  of 
appendix  B.  (Note:  The  conductive  heat  flux  at  the  front  surface, 
needed  as  the  input  for  the  Computer  Program  B,  was  the  one  predicted 
by  Computer  Program  A  of  appendix  B. )  In  general,  the  calculated 
results  are  in  good  agreement  with  the  experimental  data.  The  differ¬ 
ence  is  often  less  than  10°C,  with  a  few  exceptions  of  25°C  as  the 
maximum  difference.  However,  the  predicted  temperatures  within  the 
polyurethane  foam  show  significant  deviations  from  the  measurements. 
The  negative  deviations,  i.e.,  temperatures  measured  lower  than 
predicted,  are  seen  when  the  exposed  time  is  small  and/or 
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FIGURE  11(b)  DEPTH  FROM  FRONT  SURFACE  (cm) 


FIGURE  11(c)  DEPTH  FROM  FRONT  SURFACE  (cm) 
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FIGURE  11(d)  DEPTH  FROM  FRONT  SURFACE  (cm) 


where  the  depth  is  large.  The  positive  deviations,  i.e.,  temperatures 
measured  higher  than  predicted,  are  seen  wherever  the  void  is  observed 
inside  of  the  polyurethane  foam.  This  is  again  believed  to  be  caused 
by  a  convective  flow  of  hot  combustion  gas  within  the  void,  which 
leads  to  increased  heat  transfer  to  the  colder  side. 

The  weight  losses  calculated  by  the  model  are  compared  with  the 
experimental  data,  as  shown  in  figure  12.  The  figure  shows  that  the 
calculated  weight  loss  is  about  15  to  95  percent  higher  than  the 
measured  values  in  all  cases,  but  the  order  of  the  mass  loss  is  in 
agreement  with  the  experimental  data  except  for  the  control  (without  a 
fire  blocking  layer).  This  again  is  due  to  the  presence  of  the  void 
formed  inside  of  the  foam,  which  is  clearly  the  physical  state  this 
thermochemical  model  is  not  aimed  to  be  used  for. 

5.  SUMMARY 

The  work  on  the  prediction  of  thermochemical  performance  of  multi-  layered 

seat  cushion  materials  leads  to  four  conclusions: 

1.  The  concept  of  thermochemica 1  sublayers  is  applicable  to  describe  the 
complex  pyrolysis  behavior  observed  for  many  of  the  currently  used  and 
proposed  aircraft  interior  polymeric  materials. 

2.  The  thermochemical  model  of  the  past  can  be  extended  to  handle  such 
rjltilayered  systems  analytically. 

3.  The  model  predicts  reasonably  within  a  factor  of  2  both  the  weight  loss 
due  to  burning  and  the  temperature  profile  established  within  the  seat 
materials. 

4.  The  thermochemical  model,  can  be  used  with  the  minimum  number  of  input 
data  determined  by  experiments  for  the  thermochemical  performance 
prediction  of  other  multilayered  materials  under  fire  conditions.  This 
analysis  coupled  with  experiments  in  the  NBS  Smoke  Chamber  offers  a 
useful  small-scale  test  procedure  for  evaluating  candidate  blocking 
layer  materials. 
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APPENDIX  A  Derivation  of  r 
The  differential  equations  are 

k  —  +  Pcr  =  DPNB  exp  (-E/RT) 
dx2  dx 

.  1  dN 

Ndt  =  B  exP  <-E/RT> 

and  the  boundary  conditions  are 

at  x  =  0  T  =  Ts 

N  =  Ns  =  1  -  1 

FSs 

at  x  =  xL  T  =  Tc 

N  =  Nf  =  1  -  !_ 

FSC 

By  introducing  the  dimensionless  terms, 

y  =  Per  x 
k 

and  T  =  T-T0 

^s-^o 

Then  Eqs.  ( A- 1 )  and  (A-2)  become 

d2r  +  dr  =  kDPNB  exp  (-E/RT) 
dy?  dy  (Per)2  (Ts-T0) 

dN  .  dT  =  k  cPNB  exp  (-E/RT) 


T 


2 


Further  use  of  dimensionless  variables 


c(Ts-T0) 


convert  the  Eqs.  (A-7)  and  (A-8)  to 


(A-9) 


(A— 1 0) 


(A-ll) 


pp'  +  p  =  A  h  N  exp  ( -E/RT ) 
pN*  =  A  N  exp  (-E/RT) 


( A- 12 ) 
(A- 13) 


where  a  prime  denotes  d/dr. 


The  boundary  conditions  are  accordingly 


r  =  1 ,  p  ■  ps,  N  -  Ns 
T  -  rc»  P  =  Pc »  N  =  Nc 


(A-14) 

(A-15) 


The  Eq.  (A-l 3)  can  be  combined  to  Eq.  (A-l 2)  to  yield 


pp'  +  p  -  pN'h  =  0 


(A-16) 


Assuming  p  /  0,  then  Eq.  (A-16)  becomes  essentially 


+  1  -  dN  *  o 

dT  dr 


(A-l 7) 


If  Eq.  (A-l 7)  is  integrated  using  the  boundary  conditions  of  Eqs.  (A-14)  and 
(A-15),  then 


N  =  (p  +  T-Ps-l+h  Ns)/h 
or  N  =  (p  +  r  -  pc  -  rc  +  h  Nc)/h 


(A-l 8) 
(A-l 9 ) 


Then  substituting  Eq.  (A-19)  into  Eq.  (A-12)  gives 


pp '  +  p  =  A(p  +  r  -  ps  -  1  +  h  Ns)  exp  (-E/RT)  (A-20) 

This  is  the  differential  equation  to  be  solved  with  the  boundary  equations 
(A-l 4)  and  (A-15). 


If  the  T  is  dose  to  Ts,  then  as  an  approximation,  Eq.  (A-20)  can  be  written 
as 


PP1  + 

P  =  (p  +  r  -  ps  -  1  +  h  Ns)  A  exp  {  -  0  [1  +  X  (1-r)]  } 

(A-21) 

where 

0  =  E 

RTS’ 

(A-22) 

and 

X  =  Ts-T0 

TS 

(A-23) 

Let 

-X-  =  A  exp  (-0) 

(A— 24) 

and 

g  =  -ps  -  1  +  hNs 

(A-25) 

then. 

the  Eq.  (A-2I)  is 

pp1  +  p  =  (p  +  t  +  g)  Xexp  [-  0  X(l-r))  (A-26) 


Further  defining 

n  =  1  -  t 

e 

£  =  p  +  r 

X  =  e  -X. 

where  £  =  1 

ex 


(A-27) 

(A-28) 

(A-29) 

(A-30) 


A3 


enables  rewriting  Eq.  (A-26)  as 

-  U  -  1  +erj)  *  (  4  +  g)  A  exp  (-77) 


(A-31) 


_!L  -  ML  =  xe-7?  dr, 

i  *  9  4  ♦  9 


(A-32) 


The  boundary  conditions  for  this  equation  are 


r,*0  4  *  4S  =  ps  ♦  1 

*7  *  *?c  4  *  4  c  *  Pc  +  Tc 


(A-33) 

(A-34) 


Integrating  now  with  respect  to  4  and  H,  in  the  region  between  two  surfaces,  s 
and  c,  is 


Sc-ri  -  Sc  (?h)d«  ■  Sc  x  e*"d" 


(A-35) 


to  give  the  result  as 


(^)-l 


U,  -£«;)-  9 


(A'36) 


Then  A  Is  simply  given  by 


O'-  e-  *c)  •  }(«s  '  4C) 


+  (1  +  g)  ln^  +  g)|  (A“37) 


By  combining  Eq.  (A-29),  (A-30),  (A-22),  (A-23)  and  (A-ll),  r  can  be  shown  in 
more  familiar  terms, 


(k/Pc)  B  exp  (-E/RTS) 


(A-38) 


A4 
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Computer  Program 
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P  E  AC  IN  TPE  INGREDIENT  “ATfPIAL  PR0PEF7T  AND  PRINT  OUT 
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CO  102  I=1,1TCTAL 
NTR  =  NTREC  <  1  ) 

WRITE  <0,951)  I.ANA’U  T  < l , 1 > , A NAME T i 1  ,21 
951  FORMAT  <5» , 12 ,2> ,2A5) 

V>  1  T  f  < fc , °53 )  (ULOll  ,J>,TKH1<  1,J), 

-  KTt  1  ,  J)  ,CF  T  <  1  ,  .l>  ,KHOT  I  ]  ,  J)  ,PT  <1  ,  J>  ,ET  <1  •  J)  ,01  <  I  ,  J>  ,ME  AT^T  <  I  ,  J)  , 

J=1.\TR> 

9  53  F  OR  PA  T  I  1 7  V  ,  P7.1,  F7.1,  1*. 

F1C.5.JX,  Ffc.3,  3'*  Ffc. 3,  27,  Ell. 3,  AX,  F9.0,  F10.1,  F1C.2» 
1 C  2  CONTINUE 
C 

C  Fir  CONSTANTS  USED  IN  THIS  PROGRAM 
C 

r= i  ,o 

e=l  .9P717 
PK=«2.06 
CF tAS=C.3R2 
CPS  L  0N  =  0 .9 
$IGF'A=1«355E-12 
T  C  =  22  5, 0 
C 

C  CALCULATE  N  AND  FS  AT  TwO  END  TEMPERATURES  OF  KINETIC  RFC1“ES 
C  The  CALCULATION  SECUENTE  IS,  BY  LCGIC,  From  LOW  TO  hi  GH  TEMPERATURE 

c 

WRJTF  <6.0E5) 

955  FGPPAT  <1*M,  1C>,  ’NORMALIZED  Ffa&kENT  SIZE  *T  END  TEMPS  LT  •  ,  1», 

-  ’KINETIC  RECIPES’.  /,  11X,  5£><’:’),  /> 

DO  200  Isl.lTOTAL 

WRlTf  (6,951)  I,  ANAHETU,  1),  A  N  AM  r  T  C I  ,  2  > 

N  TR~ NINE 9(1) 

NKLO  < I.NTfi  ): 1 .0 
FKLO< I, ‘TP ) =10000. 

DO  2 C 5  JJ=1,NTR 
J=NTR’1-JJ 

CALL  FPAGPT  <PT  I  1  ,J),E It  I ,J ) , TKLOt I  , J>, TKH1 <  I, J), 

•  UEATPT ( I ,J) ,NKLO(I ,J) ,NKHI <I,J) > 

IF  (J.VF.l)  N  KL  3  < 1 , J- 1 ) =NKHI ( I , J ) 

IF  <NKLO(I,J).E0.1.0)  NKL0<l»J>=0. 9999599 
FKL0<I,J)  s  l.O/ll.O-NKLOil.JIJ 
IF  (NKHI < 1 , J) ,E a. 1. 0)  PKHII I ,J>=0. 9999999 
FKHKI.J)  =  l.OMl.O-NKHMI,  J)» 

WRITE  <E,°52)  TPL0< I,J),TKHI < 1, J),NKLO< I,J),NKM1< I ,J) 

958  FORMAT  <22*.,  F7.1,  •-*,  F7.1,  5X*  Fll.o,  2X,  •  —  •»  IX,  FU.6) 

205  CONTINUE 
2  CO  CONTINUE 

c 

C  INTT1ALL1ZE  A  FEW  VAFIABLES  FOR  COMPUTER  PLOTTING 
C 

NPAGE  =9 
NP  L  OT  =  1 


?• 

3, 

)• 

0‘ 

l* 

2« 

5* 

4* 

5‘ 
6* 
7. 
3* 
9  • 
0* 
1* 
2‘ 
3* 
4* 
f>« 
6* 
7* 
:■ 
o» 
o* 


1  PL  0 T  =  C 
C‘L  L  SC-f.PLT 
C 

C  READ  IN  TEMPERATURE  DATA,  I.E., 

C  T«-PP.  OF  TOP  SURFACE  op  ea;h  LATER  AND  REAR  SURFACE  (COLO) 

c 

1111  CONTINUE 

WRITE  <6. A  98) 

898  FORMAT  ( 1 H 1 1  25X,  •PURNINP  RATE  PREDICTION  FOR  MULT  1  -  LA  Y  ERE  0* 

-  •  AIRCRAFT  SEAT  CUSHION*  /  2EX*  634*=*>,  /> 

PE  *0  45*904*1  ND=1112>  T3ACK*(MEA01>.(1  1*  1=1*14) 

9C4  FORMAT  4F1C.0.14A5) 

■  FITE  f 6,899)  lhCADIV<I).l=l,14) 

899  FORMAT  <4  OX,  14A5,  III 

5>,  *8.  ATERIAL  *  7  A ,  ‘SURFACE  THERMAL  HEAT 

-  DENSITY  PRCt>P.  ACTIVATION  HEAT  OF  HEAT  *,/ 

-21 X ,  ‘TEMP.  CONDUCT  IV.  CAPACITY  FACTOR  ENERGY 

CECOMPOS.  RATE**/ 

-21V,  •  <  K  )  *  ,  5*,‘ICAL/CP.S.C)  (C4L/GM.C)  <CM/CMJ>»,  13Y, 
-•(CAL/MPLE  )  <CAL/CM>  (OEG  C/S£C>»  //> 

1  =  1 

1C5  CONTINUE 

READ  (9,q0E,EN0  =  llC)  IUllilHKI) 

906  FTFMAT  <12. F1J. D 
I  -  1  ‘  1 


1* 

2* 

j. 

A. 

■5* 

•6* 

7. 

•8« 
•9« 
,(/• 
11‘ 
)2‘ 
)3* 
)*• 
55* 
56  * 
57* 
58* 
59« 
VO* 
V  1  • 
V2» 
43* 
*4, 
45, 
46* 

4  7‘ 
46* 
49, 
5C« 
El* 

5  2‘ 


GO  TO  ICS 
11C  CONTINUE 
NL  A  YG= I -1 
C 

NLA  YG 1  =  RLA  V  6- 1 
DO  115  Irl.VLAYfl 
llc  TLO< I >=TH1 < 1*1 > 

TLO  <R'LAYG)  =  TPACK 

c 

r  ROW  CALCULATE  the  fragment  SIZES  AT  THE  TWO  IsOURCAKlES  OF  LaCM 
C  GENERIC  LAYER 
C 

DO  300  1  =  1  ,  F.L  A  YG 
IDENT=1D< I ) 

NT  I D  =  NTRf 0  t 1DCNT ) 

IF  (TLOU  ).GE  .TKLOUOENT  .MTIC. ) )  GO  TO  3C1 
JST  ART  < I >=NT  ID 
MOrltsl.O 
FL 0  < I >  =  10000. 

GO  TO  304 

301  CONTINUE 

00  3C2  J= 1 ,NT  ID 

IF  (TLO  <1  >  .GE.TKLO*  IOENT,  J)  .AND.UO  < 1 >,LT.T«H1 < I  DENT , J) ) 
GO  TO  305 

302  CONTINUE 
WRITE  <6,961)  J 

961  FORMAT  <  5Y, ‘FOR  •,  J2«  *TH  LAYER,  TLO  IS  OUT  OF  BOUND*) 

GO  TO  304 

305  JSTA  RT  <  I  >  =  J 

CALL  FPAGPT  CRT <1  DENT, J) »ET  < I  DEN  T , J ) , T M L 0 (1 DE NT , J) , T LO< I ) , 
MEATRTIIDENT,J),NRLOtlDENT»J),NLOM>) 

F  LOI  J  )  =  1 .0/(1 .O-NLOI I ) ) 


B3 


■3. 

4* 

•5* 

•  6* 

.6* 

:  9  • 
,G* 
»1* 
.2. 
>3« 
>4* 

■  5* 

i  6  * 

7. 

R. 

9* 

■o. 

'I* 

2* 

!• 

'S* 
'6* 
r  7  » 

n, 

(0* 

•i* 

>?• 

>3* 

i4* 

L  6* 
>7* 

Ee* 

£9* 

ec* 

9i« 

92* 

^3* 

®4« 

SS* 
5  f  « 
SI* 
9P« 
9®, 
00* 
ox* 
02* 
03* 
0  4* 
OS* 
0  6* 
C7« 
C  b* 
09* 


3C4  CONTINUE 

if  ithh  i>  .gf .  naoiiDf  a/t,ntio))  go  to  306 

J£  N  P I  1  >=NTI0 
NM  I  I  ]  >  =  1 . 0 
Fhl II >=10000. 

GO  TO  300 

306  CONTINUE 

CO  3C7  J=1.*T1D 

IF  UHI 1 1  7.SF  .TKl  GUDCUT  «J).AND.TI<I  I  I  l.LT.TKHlI  in£M,J)  ) 

-  C-  0  TO  30b 

307  CONTINUE 
VR1TC  <F,9fc2>  1 

962  FORMAT  (5X,  « FOK  *,  1?,  *Ih  INTER,  TMl  IS  OUT  OF  ROUND*! 

3  CP  JFNOIllrO 

CALL  rRAGF.T  t&TI  IDEST,jj  ,f  TM  DE>-T.  J)  ,T*LOI  IDENT.U)  ,ll«ll  l>  , 

FE  A TNT  IT  CENT ,J> ,\KLOI ICENT.J) ,Nhl tl>> 
rHlI I )r l.o/ll . C  *  N  h I  1 1) ) 

3  C  0  CONTINUE 
C 

r  RCrCFINE  AND  ACNURhER  THE  SU°D1VJDED  hULT lLAYCRS,  BASED  ON  THE 
C  TEHFEF  ATuRf  FANGS  3  CF  T«E R ROC HE R 1 C AL  FEGI"ES  AND  ALSO  ON  ThF  Of  VC  R JC 
C  DIVISIONS 

C  SftUCNCC  Or  1J  IS  FPDv  FRONT  HOTTEST)  LAYER  TO  FEAR  iCOLDEST)  LATER 
C 

IJ  =  C 

DO  310  l  =  1 *NL  AYG 
IDE  NTsl Cl  1  ) 

UENDlsJENCII) 

JST  41  s  JST  A  A  T  I  1  ) 

T  J= 1 J*1 

Tl I J) =TH1 (I) 

ANA-r I1J,1 >:A\AHET I IDtNT.l) 

ANA Kl«IJ»^)san*MI T  1 1 DENT,  2) 
rn  J)r0TllDENT«v)EN01) 

C?l I J ) =CpT IITENT  , JEND I  ) 

RMO  1 1 J)=RHPT I  I  DC  NT, JEND1) 

£I1J)  =  ETIIMFT»JFND1> 

Ell J)=OTI IDENT.JFNDI) 

KII JTrKTIICEA  T.JENDI) 

HEATRA II J>  =  HE  ATRT I 1DENT ,UEND1) 

NS<IJ>=NH1  II  > 

FSSllJIaFHHl) 

NCI  IJ)=M<LCI  ICfN’T,  JENDl  > 

FSC II  JI :FKLOt I  DENT, JENDl > 

IF  IJCNOI .EO. JST4I >  NCUJ)=NLOU> 

IF  IUENDI.E8.JSTAI)  F S C I I 0) =  F L 01 1 > 

IF  IUENDI.EO.JSTM)  GO  TO  310 
JENO l 1= Jf NO l ♦ 1 

00  311  J=JEN0I1 , JST4I 
10=1 J*1 

TUJ)=TKHI  IIOLNT.U) 

ANAdEII J.1>=ANAH[TII0E\T.1) 

ANAFE  <1 J,2)  =  AUARET 1 1 DENT, 2) 

0  1 1 J  >  =C  T IIOENT.J) 

CPI10)=CPT IIOENT.J) 

R  HO  I  1 J> =RHOT I IDENT »J) 
fiIJ>=CTI10EN'T,  J) 


B4 


1C 
1 1 
12 
13 
M 

15 

16 
17 
16 
IS 
2  C 
21 
22 
23 

2<t 

25 

26 
27 
26 

29 

30 

31 

32 

33 
xk 

35 

36 

37 
3P 
3° 
AC 
A  1 
A2 
A3 
4A 
A5 
Afc 
A7 
A6 
AS 

50 

51 

52 

53 
5A 

55 

56 
5  7 
56 
5«> 
60 
61 
62 
63 
6A 
65’ 
66 


6«IJ)  =  r-l«10CKT,JJ 
X  U  J>=n  « III  NT,  J) 

NEATRAI IJ)=HF»TA1 UPENT.J) 

MSI  1 J >  =  NKHI 1 1 DE  NT  *  J) 
f  SS«  I  J»=FHH1  IIDEMT  ,J> 

NCUJ)  =  |.KLOODCNT  ,J> 
f  SCI  IJ>:iKLOI  1DENT«J> 

IF  IJ.EC.JSTAI)  MCI  I J) -MLni T > 

IF  (J.fG.USTAl)  FSCI1 J)=FLOO> 

311  CONTINUE 

310  COM  T 1 MJE 
ML  A Y  E  R  =  1 J 
T ( 1 JO >=TbAC* 

L  6  IT  C  16.500  «l,A\AAf«I,l>tA\A6.l  U,2>»  Tl])i  Kill,  CPU).  RHOU). 

-  611).  Ell),  Oil),  lirslAUl),  I=1.MLAYFR) 

OC  format  15  X ,  1  2, 2X,  245.  ffc.l.  1  x ,  F10.5,  AX,  F6.3,  5/,  F6.3, 

E1A.3,  3V,  FA.t,  AX.F7.C.FS.3  ) 

WF1TC  16,907)  TbACK 

9  C  7  FORMAT  l°>.  •  K  E  A  F  SURF*,  F6.C»  //  > 

WRITE  lb.  910) 

510  FORMAT  15/,  •PRE01CTF0  BURNING  RATE*  ///,20X,  *6URf.  R»TC* 

-  •  SURF  A  CF  TERR.  HEAT* 

-  •  FLUX  1M  HEAT  FLUX  OUT  HEAT  C07SUHCC  F  S- TOP  SURFACE* 

•  FS-66TTC6*  /  2 1 X ,  *ICF/SEC>*.  7X ,  *1K)*,  7X, 

-  ?!  • ICAL /C62SEC  )  * .  3'  ).  //> 

X  L  A  Y 1  =  NLAYTh  -  1 

ICALC  =  0 

c 

C  CtLCULATE  the  “urviw  RATES  FOR  THE  FIRST  T1*E  ASSUMING  ThAT  THE  BURN 

C  PATE  ARE  THE  SAME  FOR  EVERY  LAYFR 

C 

DO  333  1=1,  MLAYER 
T  I  0  1  I )  =  Tin  -  TO 
Hill  =  Dili  /  (crt l  )  .  1T1D-TO) 

ALPHA ( I )  r  Kill  /  (RHOI 1 )*CPI 1 >  > 

CHI  <1 )  =  1T1I >-TC>  /T  II  ) 

THTTA1  I)  =  f I  I  )/lR*Tl 1)  > 

XIS  II  )  =  -HI  I  )  /FSSll » 

XICII)  =  -HU  >  /ASCII  ) 

TAUC  =  I T  (I  ♦  1  >  -  TC)  /ITU)  -  TO 

IF  011*1). IF. 3CC.»  GO  TO  337 

ETA  C 1  I  >  =  THETA1I)  •  CHIU)  .  U.o  -  TAUC) 

E’.C  1 1  >  =  E  XP1-ETAC  II) > 

hxics  =  ihi  i  )*xic  m>  /  (Hin  ♦  nsu)  > 

LAM  BDA  I  1  >  =  IXIStD-UCtl)  ♦  ll.C*HU  )  )  •  A  LOG  I  HXl  CS  )  I  /  1 1 .  O-EMC 11  )  ) 

PS 0  =  KlI)  •  FID  «  FXPI-THETAII))  /  1 RHO  1 1  ) » CP  1 1 ) •  LAM.BDA1I)* 
THETAU)  *  CHI1J)) 

IF  IRSO.LT.O.O  WRITE  16.95  0)  1  ,  RSO,  LA  MBDA  II  ) 

IF  IRSO.LT.O.O  RSC=-R  SO 
60  TO  500 
337  CONTINUE 

RSO  =  ALPHA  U)*bU)*EXPl-THETAU))/ITHETAU).CHIlD*Ul.  0*Pt  1  )  )• 

-  ALOGIFSSU  ) /IF  SSU) -1. 0)  )  ♦  X1S11))) 

IF  IRSQ.IT.O.C)  JUTE  I 6,° 50 )  1  .H SO »F S S U ) 

IF  IRSO.LT.O.O)  PSQ=-PSO 
5  C  0  CONTINUE 

Pull)  =  SCR  T I R  SO  ) 


85 


PS  II  )  =  -1.0  ♦  XJSI1) 

cun)  =-rs<n  •  iimn  »p« id.mioi idcpi  i> 
fci  n  =  -t sue  ♦  v i c * i > 

QOUT  I  1  >  S-PCU)*  Tioni  »  RFID  ‘RHCID  *  CPII> 

ONE 111)  =  0  H.  <  1)  -  COUT  «  1  > 

J33  CONTINUE 

c 

C  CALCULATE  PUPM'JG  RATES  S  1  P  Ui  T  A  NE  OUSL  Y  F  3  R  ALL  THE  LAYERS 

C 

3  79  CORTH'UE 

IF  lf.CALC.ro.  IOC)  GO  1  0  777 
DO  38 1  1=1.NLAYER 
RFTIII  =  F  r.  {I) 

RfRCI  1  >  =  RPH1I  *  ft HC  1 1  >  »  CPU) 
nil  =  RORCC  I  >  *  <T  I  D-T  I  1«  1)  >  ♦  D'l  >»°H01I  >*RL1«1  ) 

-  •  1 1  ./F  SSI  I  >  -  1. tf SC  I I >  > 

3F1  CONTINUE 
C  WRITE  <  fc«  *r  r  <»  > 

C  --CA  FORMAT  I  /  /  A  >.  •  *  1  •  *  i  *.  ‘PS*  .11  *  «  ♦  MS*,  ID.'f’C*.  11>«*XJC».  UN'S* 
C  -  ll>»»E»ir»»  IID'CKD  11>»»LAPBDA*) 

CO  38?  U  =  1.\L»Y1 
1  =  Mill  ♦  1  -  11 
f SI  I  I  -  0.0 
cn  spa  jzi,  'jlayer 
3PA  FSID  =  PS  I  I  >  ♦  A  »  J) 

PSID  =  /  (  T  I  C  1 1  ) »R  JPC I  1  )) 

CI1)  =  -  IISID  «  1.0  -  «!!>•<  1.0  -  1.0/FSSU)  >  > 
mil)  i  PSID  ♦  l.C 

FCIII  =  PSU1  ♦  All)/  <  T  1  0  I D  <  R (  RC I D  ) 

ricin  p  pco)  ♦  1 1 1 1  «i ) -t  o  >  n  i  ci  l  > 

LTACIl)  =  IE  I  1) /R> • ITIII.il 1*1 )  I/I Tl 11  •  T  *  1 1> 

Ef.CIl)  =  C  ypi-ETAC  f  1  »» 

GX1CS= ICI1 )♦ ‘  1CI T > >  /  CGII>  ♦  XISIl)  ) 

IF  ICXICS.LT. C.0)  WRITE  16.918)  I  •  C  X I  CS  «  C.  <  1  >  .  Y  1  C  C!)  .  « I S I  1) 

IF  (GYICS. LT.C.O)  CO  TO  777 

Ole  FORMAT  IlCX.’rviCS  I  •  »  1  2.  •  )  s*  *E  15.  6*  'XIC=*  ,1  15.6.»yiS=».C15.6* 

L  A  Ff-  0  A  C  1)  :  I/ISID-YICID  ♦  C1.C*GCI  D * A LOG  I G X 1 C S D / 1 1 . D*C NC  1 1  D 
RS3  :  ALPHA  I  1  >«BI  D‘F  XPI-THETAI  ID  /  ILAKHDAI  1  >  »THETAI  IDCUMl  D 

IF  IRSO.LT.O.O)  pso=-rso 
RB I  I )  =  SORT  I RSC  ) 

R0RCC1)  =  KB  IT)  «  R  HO  II)  *  CPI1> 

All)  =  RORC I J ) • C Tl I >-T I  1 «1 >  )  ♦  CC I ) »RHD I  I ) .RB  ll> 

•  Cl./FSSCI)  -  l./FSCCI)  ) 

C  WRITE  16.925)  I  tPSI  1)  ,»  IS  II  >.PCH  >»*IC<1>  *CI1  )  *ENCI  I  )* 

C  -  LAMPDAlD 

C  925  FORMAT  1 3X  t  12.  ei3V,E10.5)  > 

C  C I  1 1  )  =-  PSCU  *  111  D-10)*KHCl)«RHncD«CPCD 

C  OOUTID  -  -PCI  I  )»CTin-TC)  •  RPID  •RHOCl)  *  CPCD 

C  C 0 f  T  1 1  )  =  G1MD  -  OCL'TID 

3R2  CONTINUE 

C  WRITE  16.920)  II.  ANAHCH.1).  A),  APE  II  »  2  >  «  RB  I  I )  *  TCD. 

C  -  Q I  Al  I  1 )  «  OOUTID.  9NETID.  FSSII).  F  SC  I  1  >  .  I  =  1  .NL  A  YE  P  > 

C 

C  STOP  IF  CONVERGED  ENOUGH  OR  KEEP  ITERATIVE  CALCULATION 
C 

DO  386  Dl.NLAYER 

IF  lAPSlPBT  II  )/RM  D-l  .0)  .IT.  0.01)  GO  TO  386 


24* 

;C| 

?(  ♦ 
:  7* 

?e* 

5  C  • 
31* 

32* 

js« 

U| 

Jf  • 
17. 

'  w  # 

19  - 

•  0  * 
.  1* 
.  2  • 
.3* 

t  4  * 

<-6* 

*  7* 
*P* 
*9* 

s>c» 

51* 

52* 

ti* 

$<» 

*  * « 

;  fe* 
5-7* 
5°* 

to* 

(  C  • 
1 1* 
62* 

6  3* 
fck* 
6S* 
66* 
67* 
te* 
69* 

7  C  * 
71* 
7  2* 
7  3* 
7*« 
7*.* 
7t* 
77* 
7  6* 
7  =  * 
7  6  • 


’CALC  =  VCALC  ♦  1 
or,  7  C  3  79 
3P6  CONTINUE 
7  77  CC\U\UC 

kMir  <6.92P)  NCtLC 

°2?  FOk'-AT  <2>,  «*,0  OF  ITERATION  =  *.  13*  /> 

00  3PS  |  =  l,Mmfi 

onm  =  -  f  s  <  i  >  •  »t<  i  )-tc»*kp  i  n*vH7ii  >*r  i :  j 

C0u7«l)  =  -PC  !  1  >*il  (I)-TD  *  D  P  1 1  >  .RM0<],  .  PIJ> 

C  k  F  T  <  I  >  =  01UI)  -  COUT  (  1  > 
ite  coivU'iur 

*FlTt  <(*920  <1.  OK'ClI.ll*  »\A*t<I*<>*  F  t’).  7|Jt, 

UlK(l),  GOuT  <  7  )  *  3NfT<l>,  rss<l>, 

~  tC  rcovAT  cfK,i2,27.2Ai»,r  U.b,  2-<.  Fic.i.  «■«,  t».  Fir.*. 

F10.A,  O  ,  EiO.S.  (iv,  E1C.5  /> 

’50  fCRPAl  ilOX.*FSG  <**10,  •  )  =  •,  E 1 3 • £•  1C*,  t;*.  I 

ALOSS=C.C 
CO  39  3  1=1*NLAYER 

3C  3  VCrSS=k'lC9S*PMll.FKC<l>*<l./tr$U  l-l  ./F  SC  <7tl 

W  «  l  T  C  <t,r'<7)  GLOSS 

'-27  f  0  F  <*  A  T  <  /  /  V  •  -flCMT  toss  PPE0U7E3  3T  THE  ‘.P'l*.  /,}  ,*»(•;•! 
/.  HC7.riC.5v.  •  &  F  A  y  /  SC  C  .  C  k  2  OF  6<jRM\„  SuB7»0  •  ) 

C 

C  PPEllMNAM  CALCULATION  PCFOkC  PLOTTING  7  vS  X 
C 

ro  1120  1  =  1  .l.'LAYf  9 

71M7V<l>  =  «T<n-T(l*l))/riOA7|iC) 

TP  < ! , 1 )=T<  1  ) 

oo  n?i  j=i*ic 

TPtI  *  J  ♦  1 )  =  7  <7  )-Tl';TV<  1  1  *  f*  l  0  A  T  <  U  ) 

1121  C  0  V  T  1  \  U  f 
1 1  2  C  CONTINUE 

00  12 1 C  l=l,FL7Yr<< 

\SNP<  1 ,  11>=NC  <  I  ) 

00  1220  JJ  =  2  « 11 
J= 12-JJ 

iFnpii.ji.u.TfTARTiin  Nssp«i,j»  =  i.a 

IF  <  7PU  .J)  .IE  .T$TA»T<  7  H  to  To  12  21 

CALL  F  P  aGp  T  <  C  <1 ) *f  < 1) .7P  <1 . 1 1 1.1P  *  I , J ) »Hb A  7  PA  < 1 1 ,NSS»  <1 , 11) « 

-  NSSP  <  1  *  J  )  > 

1220  CONTINUE 
121C  CCfjT  I  \UC 

FHS  <  NL  A  YER  *  1 1  )  =  0. COCCI 
rc  1310  1  i  =1  .  \'L  A  YE  F 
I=NLAYCP*1-I I 
RHORP  =  RWO<  I  1  .RP  <  I  ) 

On  1320  J J  =2 •  1 1 
J=1 2- JJ 

D7  0  RK=FW0RB* <CP«1)»<TFJJ,J>-TP|1,J*1»> 

♦  D<  1  >  *< \SSP<  I  ,  J*  1> -'.SSPI  1 .  J>  >  > 

P*iS  <1  *  J>rt7DXk*FHSI  1*  J*J> 

DX0  7 1  I . JJ  =  K  <! 1 /«HSC I* U> 

1320  CONTINUE 

FPS  <1-1.1 1 >=PHS  < 1 . J) 

1310  rot  rnur 

SUP=0.C 

ro  1330  1  1=1  .  MAYER 


l 


u* 

l  =  M»Yf  P*1-I1 

P  2* 

DO  1330  JJ=1.9»2 

83* 

J  -  1  2  *  J  J 

89* 

0  V  <  I»J-2)  =  (DV'0T  <1»J)*9.*CVDT(I»J-1)*DkOT>  I  »  J*  2  >  )*TIAlT(Atl  )  /  3 

85* 

SUM  =SUK*DV  < I ,  J-2) 

66* 

X< I , J-2)=SUM 

8  7* 

1330 

CUNT1NUF 

88* 

C  CLEAR  THE  MEMORIES  CF  X»LOT  AND  TRLOT 

89* 

DO  1336  1=1,121 

9  0* 

XPL 01  II):0 . 0 

91* 

1338 

TPLOTTI  >  =  0.0 

52* 

00  134  0  1=1, SLAYER 

83* 

rr  1390  jj= 1,5 

9  4  * 

j \n  =  i  J -i >  *5  ♦  jj 

55* 

J=JJ»2-1 

96* 

>F  L  OT  < IND >  =SUM  -  x (  1 , J  > 

97* 

TF  LOT  1  HD  >  =TP  (  1  ,  J  ) 

98* 

1390 

COM  I  \UE 

99* 

I\0=!N0*1 

CO* 

>  PL OT  f 1  SC  ) = SUM 

01* 

TFL  OT  (1  r.H)  -TP  U  ,  J*2> 

02* 

IFLCT=IPLOT*l 

03* 

IF  (1PL0T.GT.5>  1PLOT  =1 

04* 

If  (lHLOT  .EG.  l.ANO. NPLOT.NE. NPAGE  >  CALL  ADVPlT 

05* 

CALL  PLFORM  ( *L  INL 1 2.5,2. 75  » 

CF* 

CALL  PLA6FL  ( «TfM°rRAluRC  PROF 1LE* ,  0. 

07* 

-  'DEPTH  FROM  SURFACE  <0*1*,  0*  • T EM PE ft  A T UP E  «K)*»  0) 

08* 

CALL  PL  SC  A  L  M 0.0,10. C> , 2, 1 OC 00 C , » 3 00 *C *  -  0 C . 01,2, 10000C> 

C9* 

CALL  PLCSI7 (0*06,0*0 

10* 

GO  TO  (13*1, 1392, 1393,1314,13951,  JPLOT 

11* 

1391 

call  ORIGIN  ( 1. C. 9 . 3> 

12* 

DC  To  1390 

3  3* 

13  92 

CALL  OPTUS  «3.0,0.G) 

19* 

DO  TO  1 39  E 

15* 

13*3 

CALL  OR  1 M  t  13.0,0.0) 

16* 

DO  TO  134*- 

17* 

1  3*4 

CALL  ORIGIN  <-6.0,-3.1> 

16* 

r-n  to  139° 

19* 

13*5 

CALL  ORIUN  «3. 0,0.0) 

2  C* 

\PLOT=NPLOT*l 

21* 

1  J9E 

CALL  PLDRAf 

22* 

CALL  P L C U * V  (YPLCT, TPLOT, 121,0, ) 

2  3* 

DC  TO  1111 

29. 

1112 

CON  T INUC 

?5» 

CALL  ENUPLT 

26* 

STOP 

27* 

F  D 

B8 


1* 

2* 

3* 

R* 

5* 

fc* 

7* 

e* 

*?« 

10* 

n* 

12* 

13* 


SUBROUTINE  FRAGVKSH*  Sf  ,  S T  I  »  S  7 F  ,  SKB  ,  S  MI  «  S NF  ) 
DOUBLE  PRECISION  S X 1 « R X2, SY 1 , S Y2» DE I 
STHI=SE/<1.9t>717*STI> 

STHF=SF/<1.9B717*STF> 

SA1=EXP(-STHF >*STF  -  E X P < -ST  Hi ) *S T  I 
SKI  =  -STHI 
SX2  =  -ST  HF 
SY1  =  PE1ISX1  ) 

SY2=DE I <SX  2 ) 

SA2=  SY2-SY1 

SNF =SNI *£XP <- ISA14SF/1 .9P717»SA2)*SB/SHR> 

RETURN 

END 
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NORMAL  I2CH  FRACPFN7  S 1U  »T  1*0  U*PS  Of  K1NE7IC  RtGl'tS 


1  COVER  f*BR 


2  VONRR 


3  IS  200 


9  PRFOK 


3  PU  FORM 
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UflfiMT  LOSS  MtlMCUO  HT  THC  FOPf  L 


I* 

2* 

3* 

4* 

5* 

£>• 

7* 

e* 

9* 

ic» 
n* 
12* 
13* 
1«* 
15* 
If* 
17* 
IP* 
lc* 
20* 
>1* 
?2* 
??» 
>4* 
25* 
>6* 
?T. 
2  H  ' 
29* 
5C* 
51* 
52* 
53* 
34* 
35* 
Sfc* 
37* 
3P  • 
39* 
40« 
»1« 
*2* 
43* 
♦  4* 
»5* 
*fc* 
47* 
4  P  • 
49. 

jo* 


C  PROGRAM  P 

C  THIS  PROGRAM  TAKES  T  HE  TEMPERATURE  AM)  C0',DUCT1VL  HEAT  FLUX  AT  THE 
C  FRONT  (HOT  3  SC P r ACE  (!.£.,  R.C.'S)  AS  THE  INPUT  DATA  AND  THEN  CALCULATES 
C  THE  n  U  k  N 1 N  G  RATE,  TEMPERATURE  PROFILE,  AND  WEIGHT  LOSS  uF  M UL T I -l A YESC 0 
C  POLYMERIC  MATERIALS  (INOLUDING  FIRE  t'L  DC  K I  \G  LAYER)  USED  FOR  AIRCRAFT 
C  SEAT  CUSHIONS. 

C  THE  DIMENSIONS  USED  IN  THIS  PROGRAM  ARE  GRA*,  GPAM-KOlE,  Ei-CELSTUS, 

C  CALORIE,  CM,  AND  SECOND. 

C 

COMMON  /TCM1 /ARAM? T ( /.2 > ,Mw(  7),ITREG(  7 > , K T ( 7 « 11  )  »CP T ( 7  , 11  )  , 

1  P  HOT  I  7,  1 1 )  ,  fT  c  7,  1 1  >  ♦  E  T  (  7,  1 1  )  ,  OT  t  7,  1 1 )  ,H  E  AT  P  T  (  7 , 1 1 )  ,  TKL  U  (  7, 1 1  )  , 

2  TKHl  4  7,11  >,TS(?(  1  >,N'(;01)  ,FS  (2C1  »  ,OS<  2C1),RP<201>.:'<201)  . 

3  NKLC(7,U  ),*  KHI  <  /,1  I  )  ,f  KLO(  7,11  >  ,FK)IJ  (7,1  1)  ,m(e>«THICK  IP)  , 

4  RHCI2C1) ,LAmPEA*R,I tluI.TC.Jl 
REAL  RT,  N,  NRH1,  N«LU,  LAMBDA, Mw 
DIMENSION  HE  ADlMlt),WLOSSI7) 

C 

C  REAP  in  THT  1NGPEDIENT  MATERIAL  PROPERTY  A  \  0  PRINT  OUT 
C 

RE  A  0  <  5  »e  0 1 )  1TOTAL 
9  Cl  f  OR  R'  AT  (12) 

PC  ICC  11=1, ITOTAL 

READ  (5,902)  1  .  AN  AHE  T  4  I  ,  1  )  ,  A  NAME  T  (I  ,2)  ,MW  ( I  )  ,  NTRE  ?’(  1 ) 

°  C2  FCRRATI I2,2A5,riC.2,12) 

NTK  =  NTR  CGI  I) 

REAP  <5,SC3MTRLOiI.J),TRH1(I,J),KT(I,J),CPT( I , J) ,PhOT ( 1  , J > , 
CT(l,J),fliI,J),tT(!,J>,HfATRT(l»J),J=l»NTR) 

9C3  FORMAT  (2FS.1,  7E10.3) 

100  continue 

WRITE  (f , 94  f ) 

94b  FORMAT  1 1 H 1  ,  DE.X,  •  bUR N  I  t.G  RaU  PREDICTION  FOR  MUL  T I -L  A  Y  n>  1 0  • 

-  •  AIRCRAFT  SCAT  CUSHION*  /  2GV  ,  631'='),  / 

-  5*.  'MATERIAL  PROPERTY*/  21X,  *TrMP.  PANGE.  THERMAL*,  f<, 

-  « HE  A  T • ,  fX.'^CNSirr  PPEC'r.  ACIIPATIOA  HEAT  OF  HEAT', 

-  5X  /  3  t>X  ,  •  C  UNDUE  T  I  V  ,  C  AT  AC  IT  Y*  ,  1  3X  ,  'FACTOR  EI.CPGY', 

-  7X,  'PCCO"PUS.  RATE',  / 

-2f>,  •  ( R )  •  ,  5X,*(CAL/CM.S.C)  (CAL/G-.C)  (GM./CMJ)',  1JX, 

-•(CAL/MOLE)  (CAl/OM)  (OEG  C/SCC)'  /) 

DO  1 02  1  =  1. ITOTAL 
NTRsNTREGI 1  ) 

WRITE  (  6,  ?  M  )  1 ,A NAMfT ( 1, 1 ) ,ANAME T (1 ,2) 

5  SI  FORMAT  (5X.I2,2X,?A5) 

WRITE  (  b  ,  9  S3  )  ITKLOII  ,J)  tTKMI  (I  ,J>  , 

-  K T ( I , J> ,CPT ( I . Jl.RHOT ( 1 ,J> ,DT ( 1 , J),ET (I ,J>,DT( 1 , J> ,HEATRT( 1 . J>. 

-  J=1,NTR) 

953  FORMAT  (17X,  F7.1,  F7.1,  IX, 

FIO.5,3',  Ff.3.  3v,  Ft.. 3.  2<»  Ell. 3,  4X,  F9.0,  F10.lt  F10.2) 

102  CONTINUE 
C 

C  FIX  CONSTANTS  USED  IN  THIS  PROGRAM. 
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c 


51* 

5?* 

53* 

54* 

55* 

56* 

57* 

56* 

5?* 

60* 

61* 

62* 

63* 

64* 

65* 

66* 

6  7* 
66  * 
69* 

7  C* 
71* 
7  2* 
73* 
74* 
75* 
76* 
77* 
7B* 
79* 
60* 
81* 
62* 
P3* 

84* 

65* 
66* 
87* 
86* 
69* 
50* 
91* 
92* 
9  3* 
44* 

95* 
96* 
97* 
9  8  * 
9  9* 
00* 
01* 
02* 
03* 
04  • 
05* 
06* 


F  =  1  .0 
"sl.56717 
Rv=82.06 
CEGAS=0.382 
f F?L0N=C.9 
$1  C- WA  =  1 .355 t - 12 
T  0  =  265. 0 
C 

C  CALCULATE  N  AND  fS  AT  TWO  END  TEMPERATURES  OF  KINETIC  KESJ-ES 
C  THE  CALCULATION  SCCUENCE  IS.  67  LOGIC*  FROM  LOW  TO  HIGH  TEMPERATURr 
C 

WRITE  t6»°55> 

955  FORMAT  »1H1,  10X,  •  I.JRI1  AL  I7ED  FRAGMENT  S12E  AT  END  T  F  MRS  OF*,  1*. 

-  ‘KINETIC  REGIMES*,  /*  11X*  S6(*=*),  /) 

DC  200  I=l*ITOTAL 

WRITE  (6*951)  .,  ANAMETfl,  3),  A NAME! (I ,2) 

N  T  R  =  N  T  R EG  11 ) 

NKLP(I.NTR)=1.0 

fklo<i,ntr>=icoog. 

DO  205  J  J 5  1 , A! T R 
J=NTR*1- JJ 

CALL  FPAGWT  < HT < I , J ) , [ T « 7 , J > , IK L 0 1 J , J > , TK HI  1  I , J ) , 

-  HfAT?T«l*J),<KLOU.J),NKHHI,J)} 

IF  €  J  .  f:  F  .  1  »  NKLO«  I,J-1»  =NKH1  <1*  J) 

IF  <  t.nil(l  ,J)  .EO.  1.0)  NKlOII *J>=0. 9999699 
FKLOU»J>  5  1.0/ll.C-NKLOCI  ,  J )  ) 

IF  (NKR  1  (  I  ,  J)  ,EQ.  1.0)  NKI.IC  I  ,J>=0. 9999599 
FKHI*I,JI  r  1.0/<1. 0-NKHI »I, J>> 

WRITE  (6,958)  7 KL C I I « J > , TKH 1 ( I , J > ,NKL 01 1 , J ) , *KH K I , J » 

958  roP-AT  (22',  F7.1,  F7.1,  5*,  Fll.b,  2X ,  •--*,  lx,  F11.6) 

205  CONTINUE 
200  CONTINUE 
C 

C  INITIALL12E  A  FEW  VAFIAPLCS  FOR  COMPUTER  PLOTTING 
C 

NF  A  GE  =9 
NP  L  OT  =  1 
IPLOT=0 
CALL  PC-NPLT 
C 

C  REAP  IN  The  INPUT  Data,  I.E.,  TOTAC  NUMBER  OT  LAYERS  AND  THEIR  CRPER 
C  THICKNESS  OF  EACH  LAVER  1C*),  AND  FRONT  SURFACE  HEAT  FLUX 
C  ANP  FRONT  SURFACE  TEMPERATURE  CK) 

C 

1112  CONTINUE 

C  CLEAR  THE  MEMORIES  OF  PLOT 
00  1336  1=1,201 
X( I >=0.0 
1338  TS ( I »=0.0 

READ  (5,907,rN0  =  1110>  < HE AD  I N ( I ) » 1  =  1 , 1 6 ) 

907  FORMAT  (16A5) 

READ  (5,904)  NLAYG, (10(1), 1=1, NLAYG) 

READ  <5,905)  (THICK (I >, 1=1, NLAYG) 

READ  (5*906)  QS1.7S1 

904  FORMAT  (1615) 

905  FORMAT  (8F1C.0I 
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)  s  • 

.9* 

.0* 

11* 

>2* 

13* 

13* 

15* 

U>* 

1  7* 
18* 
>9* 

'  I  • 

!  1  * 
’2* 
’3* 
•3* 
•5* 
;6* 

’  7* 
ie* 

2  9. 

?c« 

51* 

32* 

53* 

54* 

35* 

36* 

37* 

3P* 

39* 

30* 

31* 

3  ?« 
33* 
33* 
35* 
3fc* 
37* 
3fr* 
3Q* 

5C« 

51* 

52* 

53* 

53* 

55* 

5fc* 

5  7* 
5e* 
59* 

6  0* 
61* 
62* 
63* 
63* 


9C6  FORKATI2F10.C) 

WRITE  (6, °  59 )  I  HE*  JIM  1  >,  J=1  ,lfc> 

559  FCRRAT  (lu1.2CX»lfeAb/) 

WRITE  (  6*  S  6 0  > 

960  FORMAT  1103*  «KAT£PIAL  /  T  H I  CK  ,\  E  SS  *  *  /  > 

00  25C  1=1,NLAYG 

131=1511) 

250  WRITE  (6»°6l>  ANAKETUPJ.l),  AN  APE  T  (  1  DI  ,  2  I ,  THlCK(J) 

961  FORMAT  (  2  0  v  *  2A5*  5X*  F5.1*  »  C***/) 

WRITE  ( t  ,  9  6  3 )  CS1.TS1 

963  f  CR “A  T  ( 1 C  v  ,  ‘FOUND AT  7  CONDITION*./*  20X,  »HEAT  FLUK  (C  0  NOUC  T  I  Vt )  * 

-  *•  =  *,  F1C.3.  •  CAL/CP2.SEC*,  /«  2PX,  ‘SURFACE  TEHPERATURE  s», 

-  F 6.1*  •  K»,  /> 

C 

C  THE  OI^lNSIONS  FT R  TS*  «,  OS.  N«  RP  HAVE  TO  tE  LARGE  ENOUGH  TO  COVER  THE 
C  WHOLE  TEFRERATURf  DROP  FROR.  THE  FRONT  TO  REAR  SURFACE  IS  COVERED 
C  WITH  The  TE"PERATUFE  Interval  I T I F.  TV  )  AS  SHOWN  BELOW 

c 

1 1 \ TV  =  3 • 0 
1  =  0 

T9(1)=TS1 
C5<  1>=0S1 
' ( 1 J=0. c 
STHlCKrC.O 
TwLT=0.0 
00  251  1 N  G  =  1 , 7 

251  WLOSST  INC>=0. C 
WRITE  (  6.  °  65 ) 

°65  FORMAT  U1X,*]*,  ay,  5>*  •TSil)  *  TSlI*l)*,  8X,  *0'*, 

-  5r,  *x<n  -  v(i»i>*,  5*,»osn>  gsm*i>*,  iox,*nu>  mi*i>*/» 

DO  500  11  1  =1  ,  NLAYG 
11=1*1 
101=10(111 ) 

VP  I T  E  (6,9701  III,  A, AMET( I D1 *1  >,  ANARET (101,2) 

570  FORMAT  (3V,  )*,  i/, 2*5) 

MR=NTPEG(  101  I 
*THICK=ST*'1CK*THIC<(I  I  I  I 
5C2  CONTINUE 
1  =  1*1 

TS(  1*1  )  =  TS(  1  )  -  TINTV 
503  CONTINUE 

IF  <TS(  1*1 l.LT.TO)  TS 1 1 ♦  1  >  =T  0 
IF  (TS(1),LE  ,TKLP(I0I,MR1  )  \( I )=1 .0 
IF ( TS(I >.LE .TKLO(IOI,UTR) 1  PC  TO  310 
00  307  0=1  ,N'TR 

TF  (TS( I ).6T.TFL0(I0I ,J).ANP.TS(1 ) .LE.TKHH 1D1,J>)  60  TO  30e 

307  CONTINUE 

308  01=0 

IF  (TS( 1*1 ) .L T.TKLOd 01 *J1> >  TS(l*l»=T«Lu( 101,01) 

CALI  FRAGPT  ( R T ( 1 D 1 ♦ J ) , E T ( 1 01  , 0 > , TKLO ( 1 01 , 0  ) , T S ( I > , 

1  HEATRT (101  * 0>,NKL 0(101 ,0) ,N(I) I 

if  (N(i). ro. i.o)  rsd)  =  icoco. 

IF  (N  (I  )  ,NE  ,1.0)  FS(  I  »  =  1.G/(1.0-NU»  ) 

CALL  FRAGHT  (PT ( 101 ,J ) ,FT ( 101 ,0),TXLO( 101 .0>,TS( 1*1 )  , 

1  HEATRT  (101  ,0>,NWL  ODD  1,J),V(  1*1)  ) 

IF  (N(I*1).C0.1.0)  FS(I*1)=1CCC0. 

IF  (N(I*1).FE.1.C)  FSU*1)  =  1.C/(1.0-'J(I*1>) 
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CALCULATE  The  FUUM'.r.  RATl  R  USING  THE  HEAT  FLUX  AT  1 MF  FRONT 
AND  IrtE  ASSUMES  HEAT  FLUX  (VIA  PEAK  SURFACE  TEMPERATURE) 

CALL  HURNRT 

RHO  II  )=RMI'T  LIST  ,J1) 

DTDXraS  «I > /XT  I  101, Jl) 

px  =  cts<  D-Ts  f  i«  i  > )  /nrox 

CO  TO  40b 
310  CONTINUE 

IF  »TSCI«1).LI.TKLP(1DI.N1R)>  NUMH1.0 
OFF  PM  =  KHOT  I  1  SI  .NTT  )  *CPT  I  101  «MP  )•  <TS«  I  )-T  0) 

Rr  «  I  )  =0S«  I  )  /or  NOR. 

r  =KT(I0I,'Tc>»«TSm-TSn*l))/»EEfn*ReAI)) 

RFC  I  1  )  RKHO  T ( III, MR) 

C  S  (  I«1  >rCS  I  I )  -RHOT  I  10  1,  MR  )  «  1 )  .CPT  «  IDI  «NTP  )•  ITS  A  I  l-TS  1 1*1  )  ) 

4  OS  CCNTI*ur 

/  ( 1  *  1 ) -r I I ) *C  > 

IF  I  AMS I » I 1 *1  ) /STm I  CM -1  .C) . 17 . C.C 1 )  GO  TO  4C? 

IF  < »AS IP S *  I  ♦  1  )  )  .LC  .1 .GE-C4 )  bO  TU  499 

IF  It  PS  (IS  C 1  ♦  i  > /T  0-1. 0  )  .L  T  .0.  CC3S  >  GO  TO  4^9 

IF  IM  1  «1 )  .LT  .STHI  C«.  kor.r.S  II  ♦  1  >  .C  T.C.  0)  WM1F  16.972)  I.ROII), 
TS(I),TS<I«1),0  ',><!>.  '(I«l)(i:(II.«M*n,Ullt',(l*l) 

972  F0RMATI1CX,  1',  E1C.5,  3««  FC.l,  •  FG.l,  5V,  Ft. 2,  FR.4,*-*, 

•  F  R  .  4  ,  2  *  ,  110*4,  f  1 C  •  4 ,  2  *  ,  F1D*S,F1C*S) 

IF  IM  1*1  )  .LT  .  S1MC*.  t  NP.SSI  2*  1).  G  *.  3. 0  I  GO  TO  502 
IF  IV  l  I  *1 >  .CL .57)1  Cm)  PF  A C  T:  I  S T'i  1  C m -  <  I  1  >  >  /P x 
IF  ICSII*i >.LE  .t.C >  PFACI=0SII)/«C5I1)-3S«1*1)> 

TSI1*1)=TS«I)-ITFI1>-TS»1*1))*PFACT 
CO  TO  5  04 
499  CONTINUE 

4  RITE  1  6,  9  72)  I  ,  R  I  1  J  ,  T  M  I  >  »  I  c  I  I  «  1  >  , 

1  C  '  t  <  I  1  >  ,  X  I  I  ♦  J  )  ,  JS  I  1  )  ,  0  S I  1  ♦  1  )  ,  N  C  I )  ,  N I  I  ♦  1  ) 

,,."CTt**«  MM)  C  F  AMOVE  WILL  «l  REPLACED  PT  I  HE  FIRST  OF  THE  NEXT  LAV 
12=  I 

DC  6CC  Jr  11,12 

wLI  =RP|J)  •*  HP  |,/>  .  <•.  <  J«  1  )  *r  I  j  >  t 

lr  CWLT.LT.C.D  rfR  ITr  It  ,  *•  7 1 )  J  ,  WL  T  .  WL  US  S«  1  1  I  ) 

976  9  C  R  M  A  T  Ilf*,  MIC,  12,  *>=*.I1C.5,  •  WL  OS  S  R  *  .  F 1  0 . 5) 

t  QC  WLOSS  I  11 1 ) =  WL CSSI I  I  1 ) *»l T 

T«L  T=TwLT  *UL  CCS  I  1  I  1  ) 

5  CC  CONTINUE 

W  R  1  T F  I  6 , 9  F  2  )  I  I  M  ,  WL  i  '  S  I  1  1  I  )  .  1  I  I  :  1  LAVG  ) 

°f  2  FORMAT  I1C>,  ClbS  LOSS  FOR  •,  17,  *  L  AV  LR  =• »  Ell.*,/ 

- 1  2»  x  ,  13,  •  L*Vt Rs*,f  11,4/)) 

W°ITE  11,977)  TVLT 

977  FORMAT  12V, 51  •  •  •  ),  •  TOTAL  »ASS  LOSS  =*.F11.4,  •  G**/  SEC  »C'"2  OF  • 

-  *  Pl’RF.’ JMC  SURFACE*) 

1PLCT=1PLCT»1 
I*  IIPLOT.UT.E)  I P  L  OT : 1 

IF  I  IPLCT  .EG.  1  .AI.D.NPL  01  *NF  ,MP  AGt  >  CALL  ADVPLT 
CALL  FIFOs*  I  *L  I*.L  IN*  .2,5,2.  15) 

CALL  PLAPfL  I  *T[MprRATURC  PROFILE*.  0, 

•DEPTH  FROM  SURFACE  AC*)*.  0,  'TEMPERATURE  IK)*,  0) 

CALL  PL  SC  A  l  I  I  C.  C,  1  C.  C  I  ,2,  1  GCCC  C.  •  3CS.  0,4  CC  ,C  )  ,  2,  2  00G0C  ) 

CALL  PLCS 1 7  I C .06, L . C) 

GO  TO  I  134  1  ,  1  342, 1343 , 1  344 , 134  0) .  1PLOT 


i3*i  call  nuri\  (1.c.a.j> 

GO  TO  13 AO 

13*2  TALL  CPIGI'i  43.0.C.0) 

CO  TP  1  3*  A 

13*3  CALL  (>RJGI\  <  3.  0,  0.  0) 

GO  TO  13** 

13**  CALL  CP  I  G  J  P  4-fc.C,-3.1) 

GO  TO  13*0 

13*5  CALL  rdlll',  4  3 . 0  *  C  •  C  > 
hr  l  C 1  =  NPL  ■">  T  ♦  1 
1 3*  P  CALL  PLCPA F 

CALL  PLCUAV  4X,T«. 201,0, •*•> 
L-0  TO  1112 
111C  COA  T  !  Vl/f 

CALL  fNPPLT 
STOP 

n  r 


*  Jr  PC  uTI‘.r  f  f  AC  A  Tl  S“,S1  ,?  1  1  ,STF  ,$.tiP  ,S'Jl  ,*\F  > 
OOuLLl  PFEClJllf;  SX1,5X?,5Y1,SY2,DL I 
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PPTiPBT 1) 
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APPENDIX  C  Some  Comments  on  Burning  Tests  Procedure 


The  burning  tests  of  multi-layered  seat  cushion  assemblies  were  conducted  in  an 
NBS  Smoke  Density  Chamber.  The  orginial  radiative  heat  source  is  replaced  by  a 
Mellen  furnace  Model  10,  to  provide  a  high  heat  flux,  and  further  supplemented 
with  a  cantilever  beam  coupled  to  a  transducer,  to  measure  weight  changes  of 
suspended  samples. 

In  addition  to  these,  there  were  several  modifications  on  the  procedures  of 
sample  preparation  and  it's  burning.  These  were  instated  following  the  detec¬ 
tion  of  unsatisfactory  conditions  during  the  initial  trial  runs.  These  are 
summarized  below. 


Symptom 

Dlaqnosls 

Remedy 

A)  Uneven  heating 

1  ' 

1.  Misalignment  of 

1.  Align  correctly 

(Severe  2-D  effect) 

sample  with  furnace 

2.  Use  of  stainless 

2.  Eliminate  the  frame. 

steel  fram  as  sample 

i.e.,  large  heat 

holder 

sink  item 

3.  Loss  of  heat  on  four 

3.  Use  of  light,  thick 

sides 

Insulator  on  four 

sides,  also  as 

sample  holder 

4.  Buoyancy  of  hot  gas 

4.  No  remedy  unless 

formed  in  foam  (open 

sample  positioned 

cell) 

horizontally 

B)  Irregular  surface 

High  pressure  exerted 

Cover  fabric  cut 

shape 

by  the  gas  product 

larger  than  the 

during  Intumescence 

surface  area,  folded 

(or  disengaging  the 

over  the  sides  and 

cover  sheet  In  the 

stapled 

worst  case) 

C)  Temperature  rise  at 

Single  sheet  of  radia- 

Use  double  layers 

the  sample  surface 

tion  shield  --  not 

(with  air-gap  in 

before  removing 

enough 

1 

_ ■  ■  ■  - -  ■  J 

between)  of  shields 

Cl 


Diagnosis 


Remed1 


D) 


E) 


F) 


G) 


Symptom 


Furnace  temperature 
change  on  and  after 
removing  the  shield 


Non-linear  signal 
response  of  trans¬ 
ducer 


Ir reproducible 
weight  measurement 
of  sample 


Step  changes  In 
weight  measurement 
during  the  burning 


Change  in  Characteris¬ 
tics  of  reradiatlon 
from  the  shield  and 
the  sample  surface 


Electronic  feature 
of  the  equipment  used 


Irregularly  distri¬ 
buted  loading  on  the 
cantilever  beam 


Interference  by  the 
thermocouples 
Inserted  into  the 
sample 


Adjust  the  power 
to  the  furnace 
manual'  '  to  maintain 
the  same  furnace 
temperature 


Make  and  use  a 
calibration  chart 


Hold  the  entire 
sample  by  a  single 
wire  and  load  it  at 
the  fixed  point  of 
the  beam 


Duplicate  the  run 
without  thermo¬ 
couples  purely  for 
weight  change 
measurement 


C2 


APPENDIX  D  Photographs  of  Equipments  and  Burned  Samples 


Photo 

No.  Remarks 


1.  Thermogravl metric  Analyzer  (T6) 

The  thermograv-imetric  analyzer  is  shown  in  the  right  side  of  the  pic¬ 
ture,  connected  to  a  pair  of  flowmeters  and  a  mixer,  plus  a  vacuum 
gauge  (vacuum  pump  not  shown).  The  assemblage  of  a  chart  recorder,  a 
microprocessor,  a  heater  controller,  and  a  weight  suppressor  are  shown 
in  the  left  side  on  the  table. 

2.  Differential  Scanning  Calorimeter  (DSC)  t 

the  d§C  is  shown  here  with  a  SA t  [Scanning  Auto  Zero)  on  its  top.  The 
sample  heater  unit  is  located  on  the  deck  on  the  left  side.  The  DSC  is 
not  used  extensively  in  this  work. 

3.  NBS  Smoke  Density  Chamber 

the  h&S  Smoke  Density" Chamber  is  shown  at  left  rear.  A  newly  installed 
Mellen  furnace  is  seen  through  its  window.  The  multichannel  recorder 
(Leeds  A  Northrop,  Speedomax  Model  251)  is  shown  at  the  top  of  the 
front  cabinet. 

4.  Sample  PreparJtlon 

this  picture  shows  each  seat  cushion  sample  is  prepared  for  the  burning 
tests  in  the  NBS  Smoke  Oensity  Chamber.  The  front  cover  fabric  is  cut 
larger  than  the  face  area  and  folded  over  the  sides  and  stapled.  The 
entire  sample  is  surrounded  by  four  l"-thick  insulator  panels  (Kaowood- 
3000°  Board,  by  Pyroengineering  Co.). 

5.  Sample  Arrangement  Inside  of  the  NBS  Stroke  Density  Chamber 

Eight  pairs  of  thermocouple  (Pt  vST  Pt-ldX  Rh)  wires  are  shown  here 
penetrating  into  the  seat  cushion  assemblage,  which  is  suspended  by  a 
cantilever  beam.  The  beam  is  a  part  of  a  transducer-type  weight  mea¬ 
suring  device.  The  Mellen  furnace  is  shown  at  the  lower  left  corner. 

6.  Visual  Comparison  of  the  Burned  Seat  Cushions 

Each  seat  cushion  system  was  suddenly  exposed  to  the  Mellen  furnace 
preheated  and  maintained  at  850°C  and  then  remained  so  for  10  minutes. 
This  pictures  shows  that  the  effectiveness  of  the  fire  blocking  layer 
can  be  rated,  based  on  this  test  only,  in  decreasing  order:  LS-200 
(1/2"  thick)  >  Vonar  (3/16"  thick)  >  Preox  (11  oz/yd‘)  >  No  blocking 
layer. 

7.  Comparison  of  the  Degree  of  In-Depth  Burning 

Same  as  above,  but  this  picture  shows  the  degree  of  in-depth  burning, 
with  the  samples  oriented  in  the  same  way  as  in  the  NBS  Smoke  Density 
Chamber, 

8.  Comparison  of  LS-200  and  Vonar  as  a  Fire  Blocking  Layer 

In  this  picture,  LS-200  (s  shown  more  effective  as  a  fire  blocking 
layer  than  Vonar  when  tested  with  the  same  thickness. 
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APPENDIX  E  Comparison  With  Factory  Mutual 's  Numerical  Model 


The  thermochemical  model  (TCM)  of  this  work  is  compared  with  a  numerical  model 
developed  by  Factory  Mutual  Research  Corporation.  The  numerical  approach  is 
described  in  Appendix  A,  "A  Numerical  Model  for  One-Dimensional  Heat  Conduction 
with  Pyrolysis  in  a  Slab  of  Finite  Thickness,"  by  F.  Tamanini,  in  Technical 
Report,  FMRC  Serial  No.  21011.7  dated  November,  1976,  entitled  "The  Third 
Full-Scale  Bedroom  Fire  Test  of  the  Home  Fire  Project,"  Vol.  II,  edited  by  A.T. 
Modak. 

The  computer  calculations  were  actually  performed  with  essentially  the  same 
program  as  listed  in  the  Report,  in  order  to  enhance  the  familiarity  with  this 
numerical  model  (NM).  The  boundary  value  read-in  instructions  were  rewritten 
to  suit  this  comparison  purpose,  as  suggested  by  the  author. 

The  numerical  model  (NM)  presented  in  the  Report  is  for  a  single-layered 
material:  however,  the  extension  to  a  multi-layered  system  retaining  its 
features  seems  possible  without  difficulty.  The  change  of  total  weight  due  to 
pyrolysis  is  not  included  in  the  original  program,  but  also  can  be  added  with 
ease.  Since  such  extended/modified  computer  program  is  unavailable,  the 
comparison  of  the  result  of  the  numerical  model  with  the  experimental  data  of 
this  work  is  not  attempted. 


El 


A)  Model 


TCM 


2 

k  iLI  +  per  ~  =  DPNB  exp  (-E/RT) 
dx^  dx 

S  =  -NB  exp  (-E/RT ) 


dT  a 

^  (MgV 


3x 


apshs 

S't 


-Pa  B  exp  (-E/RT) 


For  NM,  Mg  is  the  mass  flux  of  volatiles,  h  is  the  enthalpy,  while 
subscripts  s  and  9  refer  to  the  solid  matrix  and  pyrolysis  gas  respectiely. 
Another  subscript  a  refers  to  the  unpyrolyzed  active  material. 


TCM 

NM 

Condensed  phase  reaction 

Same 

One  -  dimensional 

Same 

Steady  state 

Unsteady  state 

Moving  coordinate  system 

Stationary  coordinate 
system 

First  order  Arrhenius  reaction  for 
thermal  degradation 

Same 

Analytical  solution  available 

Numerical  solution 
technique  presented 

B)  Dimensional  Analysis 


The  dimensional  analyses  are  performed  for  these  two  models  to  determine  their 
relative  strengths  and  shortcomings.  The  dimensionless  groups  needed  to 
describe  these  models,  according  to  Pi  theorem,  could  be  the  following: 


Dimensionless  Group 

JPL's  TCM 

FMC's  NM 

Temperature 

T-  T0 

T  -  T0 

Ts  "  To 

Ts  "  T0 

Distance  Coordinate 

Prcx 

MgCX 

k 

k 

Heat  of  Pyrolysis 

D 

Q 

c(Ts  -  TJ 

c(Ts  -  T0) 

Arrhenius  Group 

E 

E 

(or  Activation  Temperature) 

~wr 

mr 

Preexponential  Factor 

Bx 

Bt 

~ 

Fourier  Number 

- 

kt 

Pcx* 

The  obvious  difference  is  that  Fourier  number,  which  compares  a  characteristic 
length  dimension  with  an  approximate  temperature-wave  penetration  depth  for  a 
given  time  t,  is  Included  in  the  NM  while  It  is  absent  in  TCM.  Other  than  the 
capability  of  the  NM  to  study  unsteady-state  behavior,  essentially  the  same 
dimensionless  groups  are  required  for  the  analysis  of  the  two  models. 


E3 


c 


1.  For  one  solid 


r 


1.  For  two  solids 

(active  material  and 
char)  and  one  Cp 
for  gas 


2.  Stepwise  variation 
with  temperature 

2.  Linear  variation 
with  temperature 

Stepwise  variation 

The  initial  and 

with  temperature 

final  densities  (char 

p 

fraction  increases 
linearly  from  0  to 
100%  of  apparent 
density) 

B.  E 

Change  often  (as 

One  set  of  constant 

evidenced  by  TGA) 

values  for  the  whole 

depending  on  tempera¬ 
ture  range 

temperature  range 

Heat  sink  term 

D:  heat  of  degradation 

Q:  heat  of  pyrolysis 

stepwise  variation 

fixed  at  a  reference 

with  temperature  (to 
be  obtained  from  OSC) 

temperature 

Boundary 

1.  Temperature  B.C.'s  at 

Either  temperature 

both  boundaries  or 

or  heat  flux  B.C.  at 
each  of  two  bounda- 

Conditions 

2.  Temperature  and  heat 
flux  B.C.'s  at  one 

ries 

boundary 

E4 


D)  Calculation  &  Results 


TCM 

NM 

Calculation 

Segment-wise  application 
of  analytical  solution, 
progressively  advancing 
from  one  boundary  to  the 
other 

Finite  difference 
representations  of  the 
partial  differential 
equation  for  the  entire 
region  and  solving  the 
resultant  tridiagonal 
matrix 

Results 

' 

surface  regression  rate, 
temperature  profile  and 
total  mass  loss 

temperature  profile, 
density  profile,  and 
mass  flux  profile 

E)  Advantages  &  Disadvantages 

TCM 

NM 

Advantage 

Parametric  analysis  is 
readily  accessible  - 
essential  for  basics 
understanding  and  estab¬ 
lishing  design  criteria 

Time-dependent  pyrolysis 
behavior  can  be  pre¬ 
dicted  (This  may  be 
useful  for  assuming 
and  confirming  time- 
dependent  boundary 
conditions  more  realis- 
tical ly. ) 

Disadvantage 

Steady  state  solution 
is  not  adequate  for  study 
of  transient  period 

Parametric  analysis  may 
be  not  as  clear  as  TCM 
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